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PREFACE 


This publication contains the proceedings of a workshop held at NASA Langley 
Research Center, January 21-23, 1982, on technology for controlling large structures 
currently being conceived for space applications. This workshop was a follow-up to 
a workshop held at Langley in October 1980. (NASA CP-2187, Structural Dynamics and 
Control of Large Space Structures (1981) contains the proceedings of that workshop.) 

The workshop focused on progress in NASA Langley’s basic research program in 
control of large space structures and heard reports on in-house efforts, university 
grants, and industry contracts, as well as on some efforts not directly supported 
by the NASA Langley Research Center program. This document contains copies of most 
of the visual material presented by each participant, together with as much descrip- 
tive material as was provided to the compiler. 

Use of trade names or names of manufacturers in this report does not constitute 
an official endorsement of such products or manufacturers, either expressed or 
implied, by the National Aeronautics and Space Administration. 


M. Larry Brumfield 
Langley Research Center 
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ACTIVE DAMPING OF A FLEXIBLE BEAM 


Garnett C. Horner 
NASA Langley Research Center 
Hampton, Virginia 



One objective in the flexible-beam research program is to develop an algorithm 
that will determine actuator and sensor locations. This is necessary because large 
space structures will have many locations where actuators can be placed. This 
research seeks to determine the "best” or optimum locations. In addition, the "best" 
locations are determined while certain constraints are satisfied which guarantee 
that mission performance requirements are achieved. 

The approach adopted in this research is to consider actuators and sensors to 
be collocated so as to produce an equivalent viscous damper. Ultimately, the 
experimental results of measuring the log decrement during free decay will correlate 
with the analytical predictions. (See fig. 1.) 


OBJECTIVE: TO DEMONSTRATE ACTIVE VISCOUS DAMPING ON A FLEXIBLE BEAM 
APPROACH: 

0 TO USE COLLOCATED ACTUATOR/SENSOR TO PRODUCE VISCOUS FORCE 
0 MEASURE DAMPING DURING FREE-DECAY 


Figure 1.- Active damping of a flexible beam. 


2 








Figure 4 shows the locations of the actuators and sensors along the beam axis. 
These locations were determined by the optimization algorithm so that the first five 
bending modes are controlled. 
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Figure 4.- Actuator and sensor locations 
on flexible beam. 
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Figure 5 shows the results of active damping as applied to the flexible beam. 

Four different sets of damper gains are shown^in the right-hand column. The vertical 

bar is proportional to the damping rate of the damper at the location shown on the 
beam. The experimental plots at the left of the figure were obtained from the first 
bending mode of the beam. With the damping rate of each damper set to zero, the beam 

was vibrated in the first mode. At the time that free decay starts, the damping rate 

of each damper is set to the desired value. .'The comparison of the percent of modal 
damping is made using two experimentally based methods. These values should be com- 
pared to the analytically determined, value. The column /entitled "graph" used a 
graphical technique to measure the amplitude of a:djacent peaks. This data is used to 
calculate the log decrement, which in turn is used to calculate percent of modal 
damping. The column entitled "ITD" contains the results of using the Ibrahim time 
domain method. There seems to be reasonable agreement between the two experimental 
methods until large modal damping is achieved. In this region small measurement 
errors can cause large modal damping errors. 


PERCENT MODAL DAMPING BEAM 



n 
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Figure 8 summarizes the accomplishments and needs of this research. 


0 ACTIVE VISCOUS DAMPING HAS BEEN DEMONSTRATED 
0 DIFFICULT TO MEASURE HIGH MODAL DAMPING 
0 ACTUATORS REQUIRE COMPENSATION 

0 DIGITAL CONTROL (VERSUS ANALOG) MAY REQUIRE MORE ANALYSIS 


Figure 8.- Summary. 
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DECOUPLING AND OBSERVATION THEORY APPLIED TO 


CONTROL OF A LONG FLEXIBLE BEAM 
IN ORBIT 


Harold A. Hamer 
NASA Langley Research Center 
Hampton, Virginia 


INTRODUCTION 


Decoupling theory is a convenient tool for devising control laws for 
structures with a large number of state variables because it allows independent 
control of each state. Complete decoupled control requires that the number of 
control actuators equal the number of modes in the system, which is a basic 
limitation in applying decoupling theory to the control of large space structures. 
Complete decoupled control is usually not achievable in practical application 
because a large space structure may have an infinite number of flexible modes; 
hence, procedures must be developed which maintain control of the structure with 
a small number of control actuators. Reduced-order systems must be utilized 
wherein only a few modes are included in the math model of the structure when 
calculating the gains for the feedback control law. In addition, some of the 
modes in the math model itself may be exempt from the control law if the number 
of actuators selected is less then the number of modeled modes. In both cases 
the control system must be designed to avoid serious problems associated with 
observation and control spillover effects caused by residual modes, which could 
result in poor performance or an unstable system. 

The present analysis presents techniques which use decoupling theory and 
state-variable feedback to control the pitch attitude and the flexible-mode 
amplitudes of a long, thin beam. An observer based on the steady state Kalman 
filter has been incorporated into the control -design procedure in order to 
estimate the values of the modal-state variables required for the feedback 
control law. 
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EQUATIONS OF MOTION 


Figure 1 shows the linearized equations o£ motion used for the decoupled- 
control, analysis of a 450-m long, thin, flexible beam in low Earth orbit. The 
equations are in modal form. The first equation represents the rigid-body 
(pitch) mode and includes the gravity-gradient effect, where is the orbital 

frequency. In the second equation n represents the number of flexible modes 
included in the math model, plus the residual modes. The damping term A 

is included inasmuch as the residual modes require a small amount of damping 
for stability. The objective is to design a control system which provides 
independent control for each of the decoupled variables. 
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BASIC EQUATIONS 


Figure 2 shows the basic equations used in the decoupled-control design. 
The equations are in state-vector form where the states x are the modal ampli- 
tudes and rates and include the residual modes. In the decoupling control law, 
u, the quantity v is the input command vector. The matrices F and G are the 
feedback and feedforward gains, respectively, which are calculated by the 
decoupling procedure. The estimator equation calculates the estimates of the 
modeled states Sc which are required by the control law. (The primes indicate 
modeled modes only.) The estimator utilizes the observation equation y and 
Kalman gains K which are precomputed by the steady state Kalman filter. The 
observation matrix senses attitude at two locations on the beam, where the 4> 
values are the corresponding slopes of the mode shapes. 
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COMPOSITE EQUATIONS 


Figure 3 shows the basic equations in composite form as a 20th order system. 
These equations are used to produce simulated time histories of system responses 
for various control input commands. The upper set of equations represents the 
system equations and includes four modeled modes and two residual modes. The 
lower set represents the estimator equations, which incorporate estimates only 
of the modeled modes. It is apparent that the control and observation spillover 
effects are caused by the matrices ^-^2x2 ‘"2x12- 

In the present analysis two torque actuators are used; each is one-sixth 
the distance from the end of the beam. One attitude sensor (e.g., star tracker) 
is at an actuator location, the other at one-third the distance from the end 
of the beam. Analyses were also performed by (1) replacing the latter sensor 
with a rate sensor, and (2) by moving this attitude sensor to the location of 
the other actuator; in both cases, however, overall performance was not as 
good as for the original setup. 


2 Control Actuators 
2 Attitude Sensors 
4 Controlled Modes 
2 Residual Modes 
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DYNAMIC CHARACTERISTICS 


The dynandcs of the system are shown in figure 4. The natural (open loop) 
frequencies and damping ratios are given in the second column. The value of 
0.001 fora), is the orbital angular velocity (orbital frequency). Small values 

O 

of damping were assumed for the flexible modes, starting at a low value of 0.005 
and increasing each successive mode by 10 percent. The fourth and fifth modes 
are taken as the residual modes. Some damping is required in these modes to 
avoid producing a system with constant oscillatory responses. 

The last two columns show values selected for the closed-loop dynamics 
for two decoupled control cases. In the FAST PITCH case the commanded pitch 
attitude is reached in about 40 seconds; in the SLOW PITCH case about 2 minutes 
are required. In the first case, the two actuators are used to decouple the pitch 
attitude and first flexible mode. The symmetric arrangement of the actuators 
produces an interaction between all four modeled modes such that a full-order F 
matrix is achieved; i.e., feedback control is available for all 4 modes in the 
math model. This condition exists only because the absolute values of the 
control-influence coefficients are the same in both columns of the B matrix. 

For other control arrangements, techniques have been developed in which the 
control-influence coefficients and/or the feedback gains are adjusted to produce 
simplified procedures for achieving overall control of the system. The current 
analysis also included model errors of up to ± 15 percent in the control- 
influence matrix, with no apparent detrimental effect on the overall system 
performance . 

In the SLOW PITCH case, the same two actuator locations were employed; 
however, the control -influence coefficients were slightly changed so that the 
decoupling control law provided control for two modes only. In this case it 
was necessary to perform two separate decoupling calculations: (1) the pitch 

and first flexible modes were decoupled, and (2) the second and third flexible 
modes were decoupled. The feedback gains obtained from both calculations were 
then combined to provide control for all four modeled modes. 


14 



q3 3 3 3 3 


DYNAMIC CHARACTERISTICS 





CLOSED- LOOP EIGENVALUES 


The symbols in figure 5 depict the eigenvalues of the closed-loop system, 
assuming perfect knowledge of the state variables. Also shown are the loci of 
closed-loop eigenvalues for the observer as the Kalman gains are increased. 

The observer (based on the steady-state Kalman filter) was designed to have a 
certain stability margin. This was attained by adding a positive scalar term 
to the diagonal elements of the system matrix A' . The gains were changed by 
varying this scalar. (The A matrix in the composite equations shown in figure 3, 
of course, is not altered.) In calculating the Kalman gains, the objective is 
to produce an observer whose response is faster than that of the closed-loop 
system with perfect knowledge of the state vector (eigenvalue real parts more 
negative than those of the closed-loop system). As shown in figure 5, there is 
no problem in meeting this condition for the filter eigenvalue corresponding 
to the pitch mode. (Hence, estimates of the pitch attitude should be very 
accurate.) However, for the flexible modes, large Kalman gains are required to 
drive the eigenvalues past the corresponding closed- loop values. In fact, it was 
found that Kalman gains which produced eigenvalues with real parts less than 
about -0.1 resulted in poor performance; i.e., excessive control forces and/or 
excessive overshoot in the flexible-mode response. An example of this is shown 
in a subsequent figure.' 
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ZERO COMMAND-FAST PITCH 


Figure 6 is an example of an instantaneous zero command (FAST PITCH case) 
to null arbitrary initial disturbances of -0.01 in pitch and O.OlS, -0.005, and 
-0.02 in the first, second, and third flexible modes, respectively, and 0.01 in 
the residual modes (4th and 5th flexible modes) . The Kalman gains used to 
determine the estimates of the modal variables (shown in 2nd column of time 
histories) correspond to real eigenvalues = -0.1, as previously discussed. It 
is assumed that instantaneous control torques T^ and T^ are available at time = 0, 

with initial estimates of 90 percent in pitch attitude and 80 percent in the three 
modeled flexible mode amplitudes; that is, the control actuators are not turned 
on until these estimates are established by the observer. Analysis has shown that 
without the effect of the controls, these estimates are achieved in about 100 
seconds. The values of initial disturbances and initial estimates quoted here 
are used for all figures which follow, except where noted. All figures except 
figure 13 pertain to the FAST PITCH case. 

As shown in figure 6, the four modeled mode responses (first column of time 
histories) are nulled after about 40 seconds. As for the residual modes, there 
is some effect on A^ during the first few seconds; however, responses in both 

modes gradually die out due to natural damping. It should be noted that, with 
the controls operating, the observer obtains good estimates of the first three 
modes, but fails in estimating A^ due to observation spillover caused by the 

residual modes. Attempts to improve this estimate by varying the Kalman filter 
gains are shown in the next two figures. 
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FILTER GAINS REDUCED 


For the zero command in figure 7 the Kalman gains were reduced so that the 
real eigenvalues, corresponding to the flexible modes, were close to zero. 

This resulted in slowly damped oscillations in the three modeled flexible modes 
due to poor estimates in all these modes. 


Time, sec 


Time, sec 


Time, sec 


Figure 7 


FILTER GAINS INCREASED 


The Kalman gains were increased for the zero command case in figure 8 so 
that the real eigenvalues, corresponding to the flexible modes, were approxi- 
mately -0.2. Here again, estimates were poor in A2 and A^. More notable is the 

poor performance as exemplified by the increased control torques and the peak 
overshoots in the flexible mode responses, which far exceeded the initial 
disturbances. Large control spillover effect is also evident in the residual 
Ag. Attempts were made to improve the performance for this case, as shown in 

the next two figures . 



Time, sec Time, sec Time, sec 


Figure 8 
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FEEDBACK GAINS ADJUSTED 


Figure 9 is similar to figure 8, except the decoupling feedback gain matrix 
was changed by deleting the gains for the three flexible- mode amplitude displace- 
ments. Hence, the control system included feedback gains only for the pitch 
attitude and rate and the three flexible- mode amplitude rates. As shown by the 
lower control torques and peak responses, some improvement in performance was 
accomplished. Also, note the large reduction in the A^ residual response. The 

system performance, however, is still unacceptable; the following figure shows 
a further attempt to improve this performance. 





Figure 9 
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1st ORDER LAG: TAU = 5 SEC 


For figure 10, in addition to the feedback gain adjustment, a first order 
lag (time constant = 5 sec) was included in the control system. This condition 
more closely resembles practical operational procedures because some lag will 
always be present in a control system. The results show further improvement 
in performance (especially in control requirements) ; however, it appears that 
the Kalman gains must be reduced for acceptable response in the flexible modes. 



Figure 10 



NO PITCH DISTURBANCE 


Figure 11 shows a zero command case similar to that of figure 6, except 
with no initial pitch disturbance. Comparison of the two figures shows the large 
effect of pitch disturbance on nulling the system. With no pitch to consider 
(and consequently no error in the initial estimate in pitch) , figure 11 shows 
considerably lower control torques and a much better response in the third 
flexible mode. 



0 20 40 60 0 20 40 60 0 20 40 60 

Time, sec Time, sec Time, sec 


Figure 11 



NO RESIDUAL DISTURBANCES 


Figure 12 is similar to figure 6, except there are no initial disturbances 
in the residual modes. Comparison of the two figures shows that the responses 
in the modeled modes are not materially affected by the motions of the residual 
modes. Also, the control requirements are about the same in both cases, except 
for the small lingering oscillations in figure 6. 






Figure 12 
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ZERO COMMAND-SLOW PITCH 


Figure 13 is an example of a zero command for the SLOW PITCH case. The 
results are similar to the FAST PITCH case (figure 6) . Although not shown, 
about two minutes are required to null the pitch attitude. Also, the maximum 
control torque is reduced by about one -half, as are the response amplitudes 
in the residual A^ . 


Time, sec 


Time, sec 


Time, sec 


Figure 13 





ONE CONTROL INOPERATIVE 


Figure 14 is an example o£ a zero command (FAST PITCH case) where one of 
the control actuators is considered to be inoperative. The feedback gains for 
the remaining actuator were not altered. The time histories show adequate 
responses in nulling the system. In contrast to the two-actuator case (figure 6), 
0 and require about three times as much time to null (Note expanded time 
scale) . 





Figure 14 
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FEEDBACK GAINS ADJUSTED, LAG INCLUDED 


Figure 15 illustrates the same case as the previous figure, except that the 
decoupling feedback gains were changed in the same manner as in the case of 
figure 9 and a first-order lag (tau = 5 sec) was included in the control system. 
Comparison with figure 14 shows better response characteristics for 9 and A^^, 

as well as a large reduction in control requirements. Also, note that the 
observer obtains good estimates of A^ after about 30 seconds. 
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TWO ZERO-COMMAND PROCEDURE 


The following four figures demonstrate a practical procedure for nulling 
initial disturbances with two separate ze ro- commands . Figure 16 represents the 
first zero-command and differs from the one in figure 6 in that actuator lag 
(tau = 5 sec) is included and the closed-loop pitch frequency has been doubled. 
Because of the increased pitch response, all disturbances are essentially nulled 
within 30 seconds. The control actuators are then turned off (observer remains 
on) at this time as shown in the next figure. 





Figure 16 


28 


CONTROLS OFF - OBSERVER ON 


The results of turning off the controls at time 30 seconds are shown in 
figure 17. The actuators were turned off in order to avoid the oscillating 
control torques which are shown to persist over a long time period in figure 16. 
As shown in figure 17, the disturbances have not been completely nulled, but 
fairly good estimates of these disturbances are obtained after 30 more seconds. 
(Note that without the controls operating the observer is able to obtain a good 
estimate of Aj.) The next step, then, is to apply the second zero command at 

this time. The resulting responses are shown in the next figure. 








Figure 17 



SECOND ZERO COMMAND AT 60 SEC 


The modal responses in figure 18 are essentially nulled after 30 seconds. 
Here, the control actuators are turned off for the final time, again to avoid 
the oscillating control torques which would be required over a long time period. 
The final results of the two- zero command procedure are shown in the next figure 
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Figure 18 







CONTROLS OFF AT 90 SEC 


The results o£ turning off the actuators after the second zero command are 
shown in figure 19. The remaining disturbances are practically zero and will 
eventually die out through natural damping. 
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Figure 19 



PITCH COMMAND - NO DISTURBANCES 


The next two figures illustrate examples of pitch commands whereby a pitch 
attitude of 0.01 radian is commanded. For the case in figure 20, there are no 
initial dusturbances, and hence no errors in the initial estimates. The commanded 
pitch attitude is reached in about 40 seconds, with only a small coupling effect 
on the second flexible mode. 
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INITIAL PITCH DISTURBANCE - ESTIMATE AT 90 PERCENT 


For the pitch command in figure 21, there is an initial pitch disturbance 
of -0.005 radian which is known only to an accuracy of 90 percent. The results 
illustrate the large effect of the initial estimate on the three modeled flexible 
modes. Also, the control -torque requirements are substantially increased over 
those of the previous figure. Doubling the error in the initial estimate would 
double the magnitude of these effects. 




03 

^ 0 




mm 



-.01 — 


1 1 

1 

-.0101 

1 1 

1 

-.02 

1 1 

20 40 

60 

0 

20 40 

60 

0 

20 40 

Time, sec 



Time, sec 



Time, sec 




Figure 21 






CONTROL REQUIREMENTS 


Figure 22 compares the maximum peak actuator torques (absolute values) required 
for pitch and zero commands, assuming no lag in the system. The data apply to a 
four-actuator control system but are representative of any control arrangement. 

For the pitch-command data, there is an initial disturbance of -0.005 radian and 
the commanded value is 0.01 radian. The zero-command data pertain only to nulling 
initial disturbance in the flexible modes; i.e., no pitch disturbance. Except 
where noted, the peak torques occurred after the initial time. 

As would be expected, the control requirements are essentially linearly related 
to the initial estimate. Also, the pitch commands require the higher control 
torques. Further, the results show that the requirements for zero commands increase 
as the accuracy of the initial estimate increase, while the opposite is true for 
pitch commands. 



Figure 22 
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SENSOR LOCATION ERROR, 5 PERCENT TOWARDS LEFT 


Figure 23 is similar to the zero command shown in figure 6, except for sensor- 
location error; i.e., the attitude sensors have been placed at locations other 
than those (nominal) used for the observer in calculating the Kalman gains. In 
this case the sensor locations have been moved 22.5 meters (5 percent of beam 
length), both in the same direction from nominal. Comparison with figure 6 shows 
that, except for A^, this error produced negligible effects on system performance. 






Figure 23 
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SENSOR LOCATION ERROR, 5 PERCENT TOWARDS ENDS 


The case in figure 24 is the same as that for the previous figure, except 
the two sensors are moved in opposite directions, where the mode-slope differences 
from nominal (for example, sign changes) are more pronounced. Even though the 
system is eventually nulled, the performance is decreased as evidenced by the 
increased oscillations in A^ and in the control actuators. These results can 

be attributed to the poor estimates in the three modeled flexible modes. 
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Figure 24 







ACTIVE STABILIZATION OF A FLEXIBLE ANTENNA FEED TOWER 


Arthur E. Bryson, Jr. 
Stanford University 
Stanford, California 


ABSTRACT 


Active stabilization logic is synthesized to hold a feed at the focus of a 
spacecraft antenna dish. The feed support structure is modeled as a tetrahedron made 
up of flexible bars and connected to the dish by six short legs containing force 
actuators. Using the symmetry of the structure, the model can be decomposed into 
four uncoupled subsystems: (1) pitch/forward motions with four degrees of freedom 

(DOF) and two controls, (2) roll/lateral motions with four DOF and two controls, 

(3) vertical motions with three DOF and one control, and (4) yaw motion with one DOF 
and one control. This greatly simplifies the synthesis of control logic. 


INTRODUCTION 

A spacecraft consists of a massive central body with a large antenna dish at one 
end; the feed for this antenna is mounted to the dish with a flexible support struc- 
ture consisting of twelve bar-like members. (See fig. 1.) Six of the bars form a 
regular tetrahedron, with the feed at the apex. Two legs connect each of the three 
joints at the base of the tetrahedron to the antenna-dish/ spacecraft , which we shall 
approximate as an inertial frame of reference due to its large mass. The mass of the 
structure will be lumped at the four joints of the tetrahedron, and the bars will be 
approximated as springs with axial deformation only. 

The design objective is to control the four lowest frequency vibration modes 
that involve lateral motions of the feed so that they are at least 10 percent 
critically damped. 


SEPARATION INTO SYMMETRIC AND ANTISYMMETRIC MOTIONS 
Motions symmetric with respect to y-z plane involve seven degrees of freedom: 

y-j^, ^2 ~ 72 “ 73 » ^2 ” ^ 3 ’ ^ 4 ’ ^4 

Motions antisymmetric with respect to y-z plane involve five degrees of freedom: 


xi, X2 = X3, Y 2 = -73, Z2 = Z3, X4 


Three of the symmetric modes involve only vertical (z]^) motions of the apex, and 
one antisymmetric mode is symmetric about the z axis (a yaw mode) , producing zero 
motion of the apex. The remaining eight modes consist of two sets of four modes that 
have identical frequencies, but one set involves symmetric motions and the other set 
involves antisymmetric motions. 

The actuator forces can be arranged into six sets, one of which controls only 
the yaw mode, another that controls only motions, and two sets of two that con- 

trol the remaining symmetric and antis 3 rmmetric modes, respectively. Thus the sta- 
bilization problem may be reduced to two almost identical problems of controlling 
four modes with two controls. 
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EQUATIONS OF MOTION 


Let be displacement vector of 1th joint, be position vector from 

ith joint to jth joint, and k 2 , k^ be equal to EA/mL^ for base members, 

vertical members, and legs, respectively. Then 


ri = -1^2™1,2“1,2' (?1 - ^2> - k2mi,3mi^3-(ri - r3> - k2mi ^ 4S1 ^ 4 • (r i - ^4) 






*2 - -''2”l,2”l,2-<f2 - - '‘l”2,3”2.3'''2 " '3> ‘ '‘A.A.i' ‘^2 ‘ 

- k£m2^5m2^5*r2 - k£m2^6m2^6-r2 + m2,5^2,5 ™2,6^2,6 




rs = -k2mi^3ini^3- (?3 - r^) - kim 2 , 3 ^ 2 ^ 3 • (r 3 - ^ 2 ) - ^ 1 ^ 3 ^ 4 ^ 3 ^ 4 • (r 3 - 

- k^m3^2“3,7*’^3,7 ~ , 8^3 , 8 ’ ^3 “3,7^3,7 ™3,8^3,8 

?4 = -k2mi,4mi^4-(?4 - r^) - ^^^2 ^ 4m2 , 4 ’ ^^4 " ^2^ " ^1®3 , 4™3 , 4 * ^^4 " ^3> 

- k£m4^gm4^9-r4 - k£m4^10™4,10-^4 + ™4,9^4,9 ■*■ “4,10^4,10 


For nominal configuration. 


= (Ddooo) 

^ (2) (10)^ 


0.5 


, = (IKIOO) _ 

( 2 ) ( 10 )^ 


0.05 


dHl . QQ) = 2.2097 
(2)(2sl2y 
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COMPUTER CODE "TETRA" 


Calculates 12x12 K matrix, where 



Calculates 7x7 Kg matrix and 5 x 5 matrix, where 




= -Va 


Ga^a 


, A r ^2 - X3 Y2 + 73 Z2 + Z3 

2 * 2 ’ 2 ’ 

, A [ ^2 + X 3 72 - 73 Z 2 - Z 3 1 

°A " [^1» 2 ’ 2 ’ 2 ’ 




T 


Note actuator forces resolved into vertical 
and horizontal components (V^, H^) at joints 
i = 2,3,4, which makes the determination of 
GgjGA quite simple. 




10 REM TETRFl **** 1/19/82 

20 REM # FINDS STIFFNESS MATRIX GIVEN 
■•0 REM # JOINT COORDINATES .ft MEMBER 
.0 REM * STIFFNESSES FOR TETRAHEDRON 
50 REM *■ ON LEGS; THEN FINDS STIFFNESS 
60 REM W. MATRICES FOR SVMMETRIC & 

70 REM * ANT I -SVMMETRIC MOTIONS. 

80 REN 

90 DIMX<10..3>.M'.’:4.. 10.SO..K<12. 12> = FORI = 1TO10 = F0RJ=1T03 : READXa .. J> :NEXTJ;. I 
1 00 READK 1 .. K2 .. KL : DEFFNR C S > = I NT C 1 0000*S+ . 5 > / 1 0000 

1 1 0 FOR I = 1 T04 •• FOR J= 1 TO 1 0 •• F0RK= 1 T03 : M ( I .• J .■ K > =X < I . K ) -X < J .• K > : NEXTK .. J . I 
1 20 FOR I = 1 T03 : FOR J=4T06 = K < I . J > =M < 1 .. 2 . 1 > #M < 1 .. 2 .. J-3 > #K2 = NEXT J I 
1 30 FOR I = 1 T03 : FOR J=?T09 : K < I . J > =M •; 1 .• 3 . 1 > *M 1 .. 3 ,J-S > #K2 NEXT J .. I 
1 40 FOR I = 1 T03 : FOR J= 1 0TO 1 2 = K <1 . J ;■ =M <1 . 4 .. I > *M < 1 .. 4 . J-9 > H«K2 : NEXT J .• I 
1 50 FOR I =4T06 : FOR J=7T09 : K i I J > =M < 2 .. 3 .. I -3 # M < 2.3. J-6 > 1 : NEXT J . I 

1 60 FOR I =4T06 : FOR J= 1 0TO 1 2 : K ( I . J > =M c; 2 . 4 . 1 -3 if:M (2.4. J-9 > #K 1 •• NEXT J . I 
170 F0RI=7T09 = FORJ=10TO12 = KCI . J>=M<:3. 4. 1 -6 > .=f:M ( 3 . 4. J-9>#K1 NEXTJ. I 
1 80 FOR I = 1 T03 : FOR J=4T0 1 2 ■■ K ( J . I > =K < I . J > : NEXT J . I 
1 90 FOR I =4T06 : FOR J=7T0 1 2 : K < J . I =K (I . J > = NEXTJ . I 
200 FOR I =7T09 = FOR J= 1 0TO 1 2 = K ( J . I > =K < I . J > = NEXTJ . I 

2 1 0 FOR I = 1 T03 : FOR J= 1 T03 K ( I . J > =-K ( I . J+3 > -K ( I . J+6 > -K ( I . J+9 > ^ NEXTJ . I 
220 F0RI=4T06 : F0RJ=4T06 '• K( I . J:;’=-K( I . J-3>-K( I . J+3;-K< I . J+6> = NEXTJ. I 
230 FOR I =4T06 : FOR J=4T06 : K C I . J > =K C I . J > -M (2.5.1 -3 *M ( 2 . 5 . J-3 i^KL : NEXTJ . I 
240 FOR I =4T06 ; FOR J=4T06 : K ( I . J > =K ( I . J > -M ( 2 . 6 . I -3 .^M ( 2 . 6 . J-3 ) #KL : NEXT J . I 
250 FOR I =7T09 ; FOR J=7T09 K ( I . J > =-K ( I . J-6 -K ( I . J-3 -K ( I . J+3 : NEXTJ . I 
260 FOR I =7T09 : FOR J=7T09 : K ( I . J > =K ( I . J ) -M ( 3 .7.1 -6 > *M ( 3 . 7 . J-6 *KL = NEXTJ . I 
270 FOR I =7T09 FOR J=7T09 : K ( I . J -K ( I . J -M ( 3 . 8 .1-6 *M ( 3 . 8 . J-6 ifKL •• NEXTJ . I 
280 FOR I = 1 0TO 1 2 : FOR J= 1 0TO 1 2 : K ( I . J =~K ( I . J-9 > -K ( I . J-6 -K ( I . J-3 .'J NEXTJ . I 
90 FOR I = 1 0TO 1 2 : FOR J= 1 0TO 1 2 : K ( I . J > =K ( I . J -M ( 4 . 9 . 1 -9 'rm ( 4.9. J-9 > .+:KL : f-^EXT J . I 
■ ...00 FOR I = 1 0TO 1 2 : FOR J= 1 0TO 1 2 : K ( I . J =K ( I . J > -M (4.1 0 . I -9 m ( 4 . 1 0 . J-9 > ■ .NEXTJ . I 

310 PRINTTAB(10>; "STIFFNESS/MASS MATRIX. UPPER LEFT QUADRANT" 

320 FOR I = 1 T06 : PR I NTTAB ( 1 9 > : FOR J= 1 T06 : PR I NTFNR ( K ( I . J : l^EXT J : PR I NT '• JEXT I 

330 PRINTTAB(10>; "UPPER RIGHT QUADRANT" : FORI=l T06 : PRINTTABC 10> ; :F0RJ=7T012 
335 PRINTFNR(K( I . J> : NEXTJ : PRINT ■■ NEXTI 
PR I NTTAB ( 1 0 ) " LONER R I OHT QUADRANT " 

FOR I =7T0 1 2 : PR I NTTAB ( 1 0 > ; : FOR J=7T0 1 2 ^ PR I NTFNR ( K ( I . J > ^ NEXTJ ^ PR I NT : NEXT I 

REN CALCULATES ANT I -SVMMETRIC. SVMMETRIC STIFFNESS MATRICES: 

DIMT(12. 12>.TI(12. 12,V.L(12. 12::-.L1(12. 12;:> : C=. 5 : T(1 . 1 > = 1 :T(2.6;;> = 1 
T(3.7> = 1 :T(4.2> = 1 :T(4.S;;' = 1 ;T(5.3> = 1 :T(5.9.:' = 1 :T(6.4;:' = 1 ••T(6. 10;:' = 1 :T(7.2>=1 
T(7.8>=-1 :T(8.9;;> = 1 :T(8.3>=-1 :T(9.10;:' = 1 : T(9. 4.'.''=-l :T(10.5;j = l : T( 1 1 . 1 1 > = 1 
T( 12. 12> = 1 : TI ( 1 . 1 > = 1 •• TI (2. 4;:'=C : TI (2. 7>=C : TI (3. 5>=C : TI (3. 8)=-C •' TI <4. 6>=C 
TI (4. 9:.<=-C ; TI (5. 10;j = l : TI (6. 2> = 1 : TI (7. 3:> = 1 : TI (8. 4>=C :TI (8. 7>=-C : TI (9. 5>=C 
T I ( 9 . 8 ) =C •• T I ( 1 0 . 6 ) =C : T I ( 1 0 . 9 =C : T I ( 1 1 . 1 1 = 1 : T I ( 1 2 . 1 2 ;' = 1 
FOR I = 1 TO 1 2 : FOR J= 1 TO 1 2 •• F0RK= 1 TO 1 2 : L 1 ( I . J =L 1 ( I . J > +K ( I . K HO' ( K . J > : NEXTK . J . I 
FOR I = 1 TO 1 2 FOR J= 1 TO 1 2 : FORK= 1 TO 1 2 : L ( I . J > =L ( I . J > +T I ( I . K > *:L 1 ( K . J : NEXTK . J . I 
PR I NTTAB ( 1 0 ; ” ANT I -SVMMETR I C ST I FFNESS.-'MASS NATR I X : " 

FOR I = 1 T05 •• PR I NTTAB ( 1 0 : FOR J= 1 T05 : PR I NTFNR ( L ( I . J > ; : NEXTJ : PR I NT : NEXT I 
PR I NTTAB ( 1 0 > " CROSS-COUPL I NO MATR I X : " 

FOR I = 1 T05 : PR I NTTAB ( 1 0 FOR J=6T0 1 2 : PR I NTFNR ( L ( I . J > ; : NEXTJ ■■ PR I NT : NEXT I 
PR I NTTAB ( 1 0 . " SVMMETR I C ST I FFNESS/M ASS MATR I X : " 

FOR I =6T0 1 2 : PR I NTTAB ( 1 0 > : FOR J=6T0 1 2 : PR I NTFNR ( L ( I . J ; : NEXTJ : PR I NT : NEXT I 
PR I NTT AE ( 1 0 > ; " CROSS-COUPL I NG NATR I X : " 

FOR I =6T0 1 2 •• PR I NTTAB ( 1 S > . : FOR J= 1 T05 : PR I HTFf-!R ( L ( I . J ; : fH'EXT J : F'R I NT : NEXT I : END 


340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

.S0i£i 

510 
520 
30 
40 
550 
560 
570 


=•'30 REM ENTER 


V.7 COORDINATES OF JOIN 

' •-•.‘-‘.'7 O 


DATA 0 . 0 . 1 v3 . 1 65 . -5 . -2 . 8b f . c; . z> 

DAT A —4 . -4 . 6 1 88 . 0 . 4 . — 4 , 6 1 88 . U . 6 . - 1 . 1 54 
REM ENTER STIFFNESS/MASS FOR HEAVV 
DATA . 5 . . 05 . 2 . 2097 


1 THRU 10: 

0 . 5 . 7735 . 2 . -6 . - 1 . 1 547 . 0 
. 0 . 2 . 5 ■ r' i' 85 . 0 . — . 5 • (' I'* 85 . 0 
MEMBERS. LIGHT MEMBERS, ft LEGS: 


READV. 
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PRINT-OUT FROM "TETRA" 


STIFFNESSVmSS MfiTRIK.. UPPER LEFT QUfiDRFlHT 
-2.5 0 0 1.25 .7218 2.0413 

0 -2.5001 -2E-04 .7218 .4167 I. 1786 

0 -2E-04 -10.0001 2.0413 1.1786 3.3334 


1.25 

• 

7218 2. 

0413 

-68. 

1694 

-14.7 

184 -2. 

.0413 

.7218 


.4167 1 

. 1 786 

-14 

a 

18<4 

-51. 

1771 

1. 1764 

2.0413 


1. 1786 

m 

34 - 

2. 

041 

3 -1. 

1764 -21.011 

UPPER RIGHT QUflliRflNT 







1.25 - 

a 

7218 -2. 

0413 

0 

0 

0 




-.7218 


.4167 1 

. 1786 

0 

1 

. 66 

67 -2 

“•S'? 


-2.0413 


1.1 786 

3.33 

34 

0 

-2. 

357 

3.3334 


50 0 


0 12.5 

21.6 

513 

0 




0 0 

0 

21.651 

.-1 

C* 0 1* 

. 502 

1 

0 




0 0 

0 

0 0 

0 







LOWER R 

IGHT QLlflliRPHT 







-68. 1694 

14.7184 2. 

0413 


12. 

5 -21 

.6513 

0 

14.718 

4 

-51.177 

1 -1. 

1764 

- 

21. 

6513 

37.5021 0 

2.0413 


-1.1764 

-21.011 

0 

0 

0 



12.5 - 


1.6513 

0 -42 

a t* r l'* 

6 

0 

0 



-21.651 

3 

■j! 1* a 50ii* 

1 0 

0 - 

76 

. 67 

09 2 

. 357 


0 0 

0 

0 li* a 

* l"' 

1.01 

1 






fiHTI-SVMMETRIC STIFFNESS/MfiSS MPTRI^' 

-2.5 2.5 1.4435 4.0825 0 

1.25 -18.1694 -14.7184 -2.0413 12.5 

.7218 -14.7184 -51.1771 -1.1764 21.6513 

2.0413 -2.0413 -1.1764 -21.011 0 

0 25 43.3025 0 -42.6776 


CRC 

ss- 

COUPLING 

MflTRIK : 


0 

0 

0 

0 0 

0 0 


0 

0 

0 

0 0 

0 0 


0 

0 

0 

0 0 

0 0 


0 

0 

0 

0 0 

0 0 


0 

0 

0 

0 0 

0 0 


SVMMETRIC 

STIFFNESS/MflSS 

MfiTRIX: 

jZl a 

5001 - 

2E-04 

1 . 4435 

0-‘0=: O 1 — -T' Ocr"7 

-2E 

-04 

-10.0001 

4. 0825 

C- m t CL O m t'OO I CL a O*-* 1 %!• a OOOT- 

a r 

218 


.0413 

—118. 1694 

-14.7184 -2.0413 21.6513 G 

.4167 

1 

. 1786 

-14.7184 

-51.1771 -1.1764 37.5021 0 

1. 

178 

6 


-2.0413 

-1.1764 -21.011 0 0 

1. 

G&G7 - 

2. 357 

43.3025 

004-3 0 —(■’6.6709 2.357 

a 

357 


a 3J34 

0 0 0 

2.3-57 -21.011 

L*RU 

SS- 

COUPLING 

MflTRD^' 


0 

0 

0 

0 0 



0 

0 

0 

0 0 



8 

0 

0 

0 0 



0 

0 

0 

0 0 



0 

0 

B 

0 0 



0 

0 

0 

0 0 



0 

0 

0 

0 0 
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EQUATIONS OF MOTION IN MODAL FORM 

Using computer code "MODALSYS**, the symmetric and antisymmetric equations 
of motion were put into modal form. Sketches of these mode shapes are given in 
figures 2 through 4. Only four modes involve fore-aft (y^) motions of the apex 
(fig. 2). Another four modes involve only lateral (x-j^) motions of the apex (fig. 3). 
Another three modes involve only vertical (zj) motions of the apex (fig. 4). One 
mode involves no motion of the apex (fig. 5). 


SEPARATION INTO FOUR SUBSYSTEMS 

Only two linear combinations of actuator forces enter into the apex motions. 
(See first example of modal controllability matrix.) We shall call them 

ffwd* different linear combinations of actuator forces enter into the x^ apex 

motions. (See second example of modal controllability matrix.) They will be referred 

to as f and f^_ . One different linear combination of actuator forces enters 

roll flat 

into the Zt apex motions. It is called f . One different linear combination of 
J- ^ vert 

actuator forces involves no apex motion, and is called f • The equations of motion 

yaw 

for these four subsystems are given elsewhere in this paper. 


ANALYSIS OF TETRAHEDRON WITH CONSTRAINED MOTION 
Symmetric Tetrahedron 


Constraints 


Xi = = 0 , x^ = -X2, 73 = y 2 > ^3 


= z. 




System equations 


X = 

(Yl, 


X2> 

u = 

(H2, 

'^2’ 


y = 

(Yl, 

Y2> 


X = 

Fx + 

Gu - 


y = 

Hx 




where 
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Units m .-sec 


IiynaJiiic: 
- 2.500 

5 Matrix 
- . 000 

F.. is: 

+ 1.443 

+ 

OO't* 

• 000 

+ 2.357 

+ 1 . fc'fc'fc' 


o=:'7 

- . 000 

-10.000 

+ 4.082 

+ 2.357 

+ 6.666 

_ 0 oc;? 

c.. ■w'O t 

+ 

3- 333 

+ .721 

+ 2.041 


-14.718 

- 2.041 

+21.651 

+ 

.000 (-118.1694) 

+ .416 

+ 1.178 

-14.718 

-51.177 

- 1.176 

+37.502 

+ 

. 000 

+ 1.178 

+ Cj CiOO 

• 'v 

- 2.041 

- 1.176 

-21.011 

+ . 0H10 

+ 

. 000 

+ 1 . 666 

d* 1 

+43. 302 

+75.004 

+ . 000 

—76. 670 

+ 

2.357 

— 2 . 35 (' 

+ -j -I 

+ . 000 

+ 

.000 

+ . 000 

+ 2.357 


1.011 

Contro 1 

Distribution Ma; 

trix 

.. G.. 

is : 





+ 

. 000 

+ . 000 + 

. 000 







+ 

. 000 

+ . 000 + 

. 000 







- 

. 500 

+ . 000 + 

. 000 







+ 

. 866 

+ . 000 + 

. 000 







•f 

. 000 

+ 1.000 + 

. 000 







+ 

. 000 

+ . 000 + 

. 000 







+ 

. 000 

+ .000 + 1 

. 000 







Feedbaol* 

; Gain Matri 

X C.. is: 








. 000 

+ . 000 + 

. 000 + . 000 

+ 

. 000 

+ 

. 000 

•f 

. 000 


. 000 

+ . 000 + 

. 000 + . 000 

+ 

. 000 

+ 

. 0100 

+ 

. 000 

+ 

. 000 

+ . 000 + 

. 000 + . 000 

+ 

. 000 


. 000 

+ 

. 000 

Outr u t D i s tr i but i on 

Matrix.* H.< : 

is' 






-f- 

1 . 000 

+ . 000 + 

. 000 + . 000 

+ 

. 000 

+ 

. 000 

+ 

. 000 

+ 

. 000 

+ 1.000 + 

. 000 + . 000 


. 000 

+ 

. 00£l 

+ 

. 000 

Hi 

sturbance Distribution Matrix 

Gfl.. 

is : 





-h 

1 . 000 

+ . 000 









. 000 

+ 1.000 








+ 

. 000 

+ . 000 









. 000 

+ . 000 








+ 

. 000 

+ . 800 








4* 

. 000 

+ . 000 








+ 

. 000 

+ . 000 









Modal Analysis 


yi 

X2 

yi 

^2 

y4 

Z4 


Ei’senM-s. lu&s sne • 

Real I mail i nary Mode no. 

-151.8372 + .0000 12 

- 85.5752 + ,0000 10 

- 23.3724 + .0000 9 

- 21.7350 + .0000 7 

- 8 . 7456 + . 0000 4 

- 7. 4730 + , 0000 3 

1 . 8009 + . 0000 1 


E i sieri'v'ector Matr i x .< T .■ i s • 



. 00010 

- 

.0317 

+ 

. 0000 

+ 

. 1979 

-- 

. 2629 

- . 0000 

+ 1 

. 0000 

-f 

.0512 

- 

. 00O0 

- 

. 7307 

- 

. 0001 

- 

. 0008 

+ 1 . 0000 

- 

. 0000 

— 

. 8659 

+ 

. 9673 

- 

. 0275 

+ 

. 005 1‘‘ 

+ 

. 0832 

+ .0107 

+ 

.0165 


. 50001 


. 6753 

- 

. 0160 

- 

.1190 

+ 

. 8557 

+ .0064 

+ 

. 0997 

— 

.0193 

+ 

.0188 

+ 

. 9994 

- 

. 4997 

- 

. 1214 

+ .2440 

+ 

. 0534 

+ 1 

. 0000 

+ 1 

. 0000 

+ 

. 0316 

- 

. 1090 

+ 1 

. 0000 

- .0120 

+ 

. 1283 

- 

.0193 

- 

, 0376 

+ 1 

. 0000 

+ 1 

. 0000 

+ 

. 2424 

+ .2441 

- 

. 1 069 


®12 

®10 

mg 

my 

^4 

m 3 

mi 
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Of 

^ -the Ei-; 

i<eri*-.-‘ector M-atr : 

L X i £ : 





+ . 0000 

f 

. 0170 

- 

. 57 66 — > 3329 

- .0128 

+ .3329 

- 

. 0064 


- .0083 

+ 

. 00001 

4- 

.51 £19 - .3566 

+ . 0099 

+ . 2640 

- 

. 0099 


- . 0000 

- 

> 206 f' 

- 

.01-55 - .0090 

+ . 5656 

+ . 0089 

4 

. 2826 


+ .1253 


. 010010 

4- 

.0072 - .1-507 

- .6334 

- - 0690 

4 

.6331 


- .0997 

- 

. 0003 

+ 

.0631 + .6493 

- .0921 

+ .3793 

4 

. 0920 


- .0000 

+ 

. 8480 

4- 

.0182 + .0109 

+ .4138 

- .0102 

4 

. 2070 


+ .9486 

— 

. 00010 

4* 

.0313 + .1892 

+ .1014 

+ .1218 

— 

. 1014 


The Mod-: 

5.1 

Iiyn-3JTi i 

O: 

s Matrix Fq , Iru 

-•a)*F*iT.. 

is • 




-*•. 

+ 

. 00100 

4 

. 0000 + .0000 

+ . 0005 

+ . 0000 

4 

. 0000 

(-151.8) 

+ . 0000 


i.mm 

- 

. 0000 - . 0000 

+ .0005 

+ .0000 

4 

. 0000 

(-85.57) 

+ .0000 

+ 

. 01000 


- .0003 

+ .0000 

- .0002 

4 

. 0000 

(-23.37) 

+ .0000 

+ 

. 0010101 

4 

.0003 

+ . 0002 

+ . 0000 

- 

. 0000 

(-21.73) 

+ . 0000 

f 

. 01000 

4 

. 0000 + . 0000 

“8. 7466 

+ . 0000 

4 

. 0000 


+ . 00001 

+ 

. 0101010 

- 

.0002 - .0002 

- . 0000 

-7.4721 

4 

. 0000 


+ . 0000 

+ 

. 0000 

- 

. 0000 - . 0000 

- . 0002 

+ . 0000 

-1 

. 8009 


The Mod- 

5.1 

Corrtre 

) 1 l-Bh i 1 i ty M-atr i : 

X G , In’-.-"' 

.T>5|tG.. is 




- . 0000 

- 

.0128 

- 

. 0064 

q 





- .5643 


. 0099 


. 0099 






- . 0000 

+ 

. -5656 

4 

, 2826 






- . 1341 

- 

. 6334 

4 

.6331 






+ . 53y7 

- 

.0921 

4 

. 0920 






+ . 0004 

+ 

.4138 

4 

. 2070 






+ . 1 482 

f 

. 1014 

- 

. 1014 






The Modi 

5l 

Obsero-yhi li"ty Matrix 


is ■ 




+ . 01000 

- 

«031? 

4 

.0000 + .1979 

^ .2629 

- . 0000 

f 1 

. 0000 


+ .0512 

- 

- 0000 

- 

.7307 - .0001 

- . 0008 

+1 . 0000 

- 

. 0000 


The Mod. 

5l 

Distur 

-b. 

sbi lity M-atr ix 

G, , Inv 


.• i 

.s : 


+ . 00010 

+ 

.0170 



Aq 





- . 0083 

4- 

. 0000 








- . 0000 

- 

0 2067 








+ . 1253 

4- 

. 0000 








- .0997 

- 

. 0003 








- . 0000 

4* 

. 84S01 








+ .9486 

- 

. 00100 









Ant is3rmme trie Tetrahedron 

Constraints 

X3 = 73 - Z3 = -Z2, = 0 , = 0 , y^ = 0 , z^ = 0 , 

II3 = h2, V3 = -V2. = 0 

System equations 

X = (x^, X 2 , y2> Z 2 ’ ^4^ 

u = (h2, V2, h^) 
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Modal amplitudes 


^ t®i2» 

q = Fq + Gu + G. V 
^ q Aq 

y = H^q 


where 


Un i ts m .. sec 


Dynejfi i ci 

? Matrix F.* is- 




- 2.500 

+ 2.500 + 1.443 + 

4.082 

+ 

.000 

+ 1.250 

-18.169 -14.718 - 

2.041 

+ 12 

.500 

+ .721 

-14.718 -51.177 - 

1.176 

+21 

.651 

+ 2.041 

— 2 . 04 1 — 1 a 1 f' 6 —2 

1.011 

+ 

. 000 

+ . 000 

+25.0001 +43.302 + 

. 0010 

-42 

. 677 

Contro 1 

D i str ibut i on Matr i 

,* (j* 

is- 


+ . 000 

+ .0100 + .0010 




. 500 

+ . 000 + . 000 




+ . 866 

+ .0100 + .0100 




+ . 01010 

+ 1 . 0100 + . 0010 




+ . 0100 

+ .000 + 1.000 




Feed beet': 

C+ain Matrix C.. is 




+ . 000 

+ . 000 + . 000 + 

. 000 

+ 

. 000 

+ . 0010 

+ . 0010 + . 000 + 

. 000 

+ 

. 000 

+ . 800 

+ .0100 + .0100 + 

. 000 

+ 

. 000 

Outrut Distribution Matrix 

.. H.. i 

s * 


+ 1.000 

+ „ 000 + . 000 + 

. 0010 

+ 

. 000 


Disturb.ance Distribution Matrix Gfl^ is = 

-f- 1 . 000 
+ . 000 
+ . 000 
+ . 0100 
+ . 0100 


^1 

X2 

Y2 

Z2 

24 


E i senue. lues ere ■' 


R^.=, I 

Imayiner-y 

Mode no. 

- 85.5753 

+ , 00010 

11 


- 21.7350 

+ . 0000 

8 


— 1 1 * • 6 i‘‘ r* 5 

+ . 00001 

6 


- 8. 7462 

+ . 0000 

5 


- 1 . 8009 

+ . 0000 

2 


Eieenuector Matrix.- T.- is-' 



+ » 0256 - 

.2286 - .0001 - 

.2761 +1 

. 00100 

- .3579 + 

.1300 - .4998 +1 

. 0000 + 

.1188 

- .7840 - 

.0068 + .8659 + 

.0874 + 

.0165 

- - 0264 + 1 

, 0000 + . 0004 - 

.2207 + 

. 0926 

+ 1.0000 + 

.1409 +1.0000 + 

.8483 + 

.0901 

I n verse o*f' 

the E i •senuector 

M-atrix 

is : 

+ .0103 - 

.2877 - .6303 - 

.0212 + 

.4019 

- . 1 085 + 

.1234 - .0065 + 

« 9496 + 

. 0669 



- . 0000 - 

.3:333 + 

.5774 + 

. 0002 

+ 

. 3334 


- . 0949 + 

.6875 + 

.0601 - 

.1518 

+ 

.2916 


+ .9487 + 

.2254 + 

. 0 : 3 1 + 

. 1757 

+ 

. 0855 


The Modal 

Iiyna/dics 

: Matrix.. 

Ir'v'<T>sfF.*(fT.. is: 


+ 

. 0000 + 

.0000 + 

. 0000 

+ 

. 0000 

(-85.58) 

+ . 0000 ~'i 


.0000 + 

. 0000 

+ 

.0000 

(-21.73) 

+ .0000 + 

.0000 


. 0000 

+ 

. 0000 

(-17.68) 

+ .0000 + 

. 0000 + 

.0000 -£ 

! . i' 4fc'2 

+ 

. 0000 

+ . 01000 + 

.0000 + 

.0000 + 

. 0000 


. 8009 


The Moda 1 

Coritro 1 labi lity 

Matr i .'i 

c.< 

Iriu<T>#C 

J.. is : 

- .40120 - 

.0212 + 

.4019 





“ . 0673 + 

.9496 + 

. 0669 





+ . 6667 + 

.0002 + 

.3334 





- .2917 - 

.1518 + 

.2916 





- ,0855 + 

.1757 + 

.0855 





The Modal 

Obseruabi lity Matrix.. 

HsffT.. is: 


+ . 02.56 - 


.0001 - 

.2761 

+1 

. 0000 



The Modal DisturbalDi lity Ma'fcrix Inot i is - 

+ .0103 

- . 1085 

- . 0000 

- .0949 
+ . 9487 

Residues.: fl 1 1 Outputs.. Ctrl l..Theri Ctrl 2, etc. 

- .0103 + .0154 - .0000 + .0805 - .0855 

- .0005 - .2171 + .0000 + .0419 + .1757 

+ .0103 - .0153 - .0000 - .0805 + .0855 

Res i dues.: fl 1 1 Outputs.- Dist. 1.- Ther Dist. 2/ etc. 
+ .0002 + .0248 + .0000 + .0262 + .9487 

PITCH/FORWARD TRANSLATION SUBSYSTEM 

= -(1.342)^m^ + 0.1907fp 

= -(2.957)^m^ + 0.6652fp 

m_ = -(4.662)^m, - 0.9848f + 0.7914f^ 

7 7 p F 

m,- = -(9.251)Vn - 0.1320f - 0.5788f_ 

10 10 p F 

where 


^1 


”1 

-0.2629 

0.1979 

-0.0317~ 



^1 


0 

0 

0 

0 


”4 

^2 


0.0165 

0.0832 

0.0057 

0.9673 


m? 

^2 

= 

0.0997 

0.8557 

-0.1190 

-0.6753 


_”io 

"2 


0.0534 

-0.1214 

-0.5000 

0.0188 



^4 


0.1283 

1 

-0.1090 

1 




-0.1069 

0.02424 

1 

-0.0376 
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and 




ROLL/LATERAL TRANSLATION SUBSYSTEM 
= -a.342)^m„ + 0.2202f 

^ K. 

= -(2.957)^m + 0.5482f 

O Xj 

= -(4.662)^m„ + 0.9845f - 0.5925f 

o K L 

= -(9.251)^m + 0.1880f_ + 0.6184f 

1 J. K L 



1 

0.1188 

0.0165 

0.0926 

0.0901 


-0.2761 

1 

0.0874 

-0.2207 

0.8483 


-0.2286 

0.1300 

-0.0068 

1 

0.1409 


0.0256 

-0.3579 

-0.7840 

-0.0264 

1 




VERTICAL SUBSYSTEM 

= -(2.734)^ni2 + 0.6208f^ 
fflg = -(4.835)^mg + 0.8476f^ 

= -(12.322)^nij^2 + 0.0192f^ 

where 


50 



V 


~0 

0 

0 



1 

-0.7307 

0.0512 

^2 


0.0107 

-0.0275 

-0.8659 

^2 

= 

0.0064 

-0.0160 

-0.5000 

^2 


0.2440 

0.9994 

-0.0193 

^4 


-0.0120 

0.0316 

1 



0.2440 

1 

-0.0193 


and 



YAW SUBSYSTEM 

Jig = -(4.204)^Big + l.OOOfY 


where 




" 0 



-0.5000 

^2 

= 

0.8660 

^2 


0 

""4 


1.000 



— 


and 




1 

2 



«3 

= 

1 


1 

1 





= 0 
= 0 
= 0 



Figure 6 shows the combinations of controls that control only modes 1 and 4 
(and also modes 7 and 10) . Figure 7 shows the combinations of controls that con- 
trol modes 2 and 5 (and also 8 and 11) . Figure 8 shows the combinations of con- 
trols that control the vertical and yaw modes. 
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Top view 


Figure !•- Antenna feed tower 












Rear view 


Top view 


(a) Mode (roll) controls (f^) . 



Rear view Top view 


(b) Mode ra^ (lateral translation) controls (fj^) • 


Figure 7.- Roll and lateral translation controls. 




Side view 


Hear view 


(a) Vertical modes control (f ) . 

V 



Side view Rear view 


(b) Yaw mode controls (f ) . 


Figure 8,~ Vertical and yaw controls. 



ACTIVE CONTROL OF FLEXIBLE STRUCTURES WITH SEPARATED SENSORS AND ACTUATORS 



R. H. Cannon and D. Rosenthal 
Stanford University 
Stanford , California 
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KEY CONTROL PROBLEMS OF FLEXIBLE SPACECRAFT 


Very flexible 
Low damping 
Noncolocated sensor 

Kalman filter 


• Uncertain parameters 

• Nearly equal modes 

and actuator 

APPROACH TO CONTROL 

• Adaptive control 

Identify plant frequency 
Tune Kalman filter 





Initial laboratory two-disk plant. 
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(a) Naturoil frequency modeled 10 percent lower than actual 
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(b) Natural frequency modeled 20 percent higher than actual. 

Step response of closed-loop system with mismatch in actual and modeled natural 
frequency. 
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Step response of closed— loop system with self— tuning adaptive optimal controller. 
Vibration frequency in Kalman filter was incorrectly assumed to be 1.7 Hz (—25 
percent error) initially. Parameter error in model was corrected in real time by 
frequency-locked loop. 





DOCKING OF A RIGID SPACECRAFT WITH 
UNRESTRAINED ORBITING FLEXIBLE STRUCTURE 


T. R. Kane 
Stanford University 
Palo Alto, California 


Purpose : 

Establish a methodology for producing simulations of the docking of a 
rigid spacecraft with any unrestrained flexible structure. 

Method: 


Produce a detailed algorithm for the simulation of the docking of a rigid 
body with a nonuniform unrestrained beam. 

Background : 

Predocking and postdocking motions can be simulated with algorithms given 
In references 1 and 2. 

Tasks : 


(1) Determine the values acquired during docking by the variables character- 
izing the motions of the rigid body and the unrestrained beam. 

(2) Transform to the variables used In connection with the cantilever beam. 

(3) Perform simulations of the total process. 


REFERENCES 


1. Levinson, David A.; and Kane, Thomas R. : Simulation of Large Motions of 

Nonuniform Beams In Orbit. Part 1 - The Cantilever Beam. J. Astronaut. Scl., 
vol. 29, no. 3, July-September 1981, pp. 213-244. 


2 . 


Kane, Thomas R. ; and Levinson, David A.: 
Nonuniform Beams In Orbit. Part 2 - The 
vol. 29, no. 3, July-September 1981, pp. 


Simulation of Large Motions of 
Unrestrained Beam. J. Astronaut. 
245-275. 


Scl. , 
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0.0 iboo.o 2000.0 3000.0 ^000.0 

t, sec 

Libration of initially undeformed beam; imperfect docking. Velocity mismatch 
approximately 1 cm/sec; angular velocity mismatch approximately 0.01 rad /sec. 





o 



o 



t, sec 

Libration of initially deformed beam; perfect docking. 
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RESEARCH ON ELASTIC LARGE SPACE STRUCTURES 


AS "PLANTS" FOR ACTIVE CONTROL 


H. Ashley and A. von Flotow 

Stanford University 
Stanford, California 



INTRODUCTION 


• Modeling of large space structures (LSS) in terms of elastic wave propagation 

• Scale effects on structural damping 

• "Loss coefficients" of monolithic LSS 

• Wave propagation In nondlsperslve and dispersive media (1-D and 2->D) 

• Spectral separation of system response: 

x(s) = H(s)[u(s) + iC(s)] = Hjj(s)[u(s) + lC(s)] 

+ [H(s) - Hj^(s)][up(s) + iCp(s)] 

correction term 

where H (s) is a reduced-order transfer function 
R 

• Reflection of waves from boundaries 

• Modeling of discrete structures as equivalent continuous structures 

• Dynamics of networks of elastic waveguides 

• Control of systems with wave-related time delays 

•Application to a 1-D system under active control: 0.12-sec lag predicted with 

Timoshenko beam idealization and empirically determined shear rigidity 

•Significance of passive damping (ref. 1): 

1. A LSS with exactly zero damping is uncontrollable unless sensors and 

actuators are all collocated (often impractical) 

2. Even very small amounts of damping are important to practical success 

of control 

• Some approximate effects of LSS linear scale L on a typical modal damping 

ratio C* 

1. For a "monolithic" element, 5 is proportional to material damping and 

decreases with decreasing frequency (a) (i.e., with increasing L) 

2. Viscous friction dominates at joints; thus C 1/L 

3. Coulomb frection at joints and joint preload is dependent on rotational 

rate (nL)^ 

4. All sources active C between a constant and ^L ^ 
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STUDY OF INTRINSIC DAMPING IN MONOLITHIC METALLIC STRUCTURE 


• Two ”seTni-*reversible” mechanisms seem feasible for LSS: 

1. Thermal relaxation 

2. Grain boundary relaxation (can give large values of C but required 

temperatures may be too high) 

• Work in progress on thermal damping 
•Properties of thermal damping 

1. Involves coupling between mechanical and entropy waves; e^g., ^or ^ 
isotropic solid with T = T + AT and displacement = ui + vj + wk. 


k v2(at) - c ^ - V— 

p V 3t p[l - 2v] 3t 


• ^ = 0 


= E = aAT + 
3x XX 


o - v(o + a ) 

XX yy zz 


3a 

XX 

3x 


3a 

3y 


3a 

xz 

3z 


3t^ 


-2 -3 

2. Z is configuration-dependent (e.g., 10 to 10 for beams and plates, 

“7 “8 

10 to 10 for bars and rods). Composite beams are under study. 

3. The value of C depends on frequency O) and material properties. E.g., for 

a rectangular beam of depth b: 
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Metal (R.TO 


A1 and alloys 


Low-carbon steel 
Ti and alloys 
Ni and alloys 
Be 
Mg 

A1 at 1000 K 


b, cm 



£ r — 
GA^ = 


<0. ^2-S- 

X, 7e^/o~^ 

O.S'JS- /y-^ ^ 

30,7 /V' T^Af/JcA-rf^*/) 


fO S,6 3.2. /•& f O.J3 0./e 

WA VE'L £"//<r 7~^ 


Timoshenko beam waves. 










Thermal damping theory compared with recent tests on 
free-free A1 beams in vacuo. 



IDENTIFICATION AND CONTROL OF SPACECRAFT 


C . S . Greene 
M, F. Barrett 

Honeywell Systems and Research 
Minneapolis , Minnesota 



IDENTIFICATION & CONTROL OF SPACECRAFT 


* THE PROBLEM 

* CONTROL 

- CLASSICAL 

- MODERN 

* IDENTIFICATION 


4 z 

I 

I 

l 



ARM 
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CLASSICAL CONTROL 


* SCALAR FEEDBACK DESIGN 


* LEAD-LAG CONTROLLER 


K(S) 


k (s-fa) 
(s^+ 25 b+b^) 


* BANDWIDTH AROUND 1 RAD/SEC 


SHUTTLE AND PAYLOAD: FROM ALPHA TO Y-ATTITUOE 



FREQUENCY <rad/'scc> 



le m CL 


SHUTTLE AND PAYLOAD: FROH alpha to y-attitude 



MODERN CONTROL 


* START WITH FULL-STATE DESIGN 

- GOAL: MINIMIZE PAYLOAD ATTITUDE ERRORS 

- ITERATE ON CONTROL PENALTY TO ACHIEVE BW OF 5 r/s 

* DESIGN FILTER TO RECOVER LQ RESPONSE 

- USE STEIN/DOYLE ROBUSTNESS RECOVERY RESULTS 

* TEST ROBUSTNESS 
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1o0 frequency 


plant. r: Singular values of K<s>*C*Cs> q»lE8 





IDENTIFICATION 


TEST 

SIGNAL 

U 


n[uuT*^ 


MAXIMUM LIKELIHOOD ESTIMATION (MLE) 



MODEL STRUCTURE 


• State Space 

n* Modes (2x2 blocks) 


OUTPUTS 


t Frequency Domain 


(M= 


I 0 1 


jl. 


(t)+ 


[uW+w(^ 




Ci 0 ; , 




['] 


(t) 


y(s) = (GJs) lu(s)+w(s) 1 + v(s) 


r ' 

GJ[s) = V - 2 -^-^ -2 

t Parameter Vector 

“* - 2C*^co#^, b*^, c*^; i = 1, n*) 


85 







MLE IDENTIFICATION SETUP 


Residual Definition 

rj^ A y(kT) - y(kT) 

Likelihood Function (Negative Log) 

A N 

L(a) = Z L. (a) 
k=0 

A ** 1 r T 1 

= ^Z^ 5 [log DET S,^ r,^J 

WHERE 

a = (af. b^, i = 1, n) = UNKNOWN PARAMETERS 

Sj^ = E^{rj^r^) = predicted RESIDUAL COVARIANCE 


HE IDENTIFICATION SOLUTION 

# Parameter Estimate (Theoretical) 

a = ArG I L(»)| = PARAMETER ESTIMATE 


• Iterative Algorithms 

Gradient: Oj+j = “j - ''L(Oj) 

Newton-Rhapson: “j+2 “ “j ' [’^L(»j)] ^vL(“j) 

VL(a) ^ (a) 


WHERE 



IDENTIFICATION ACCURACY ISSUES 


• SYSTEMATIC ERRORS: £*(&) - a* 

- MODEL-ORDER MISMATCH 

- TEST SIGNAL AMPLITUDE AND SHAPING 

- SYSTEMATIC DISTURBANCES 
-SENSOR/ACTUATOR MODEL ERRORS 

• STOCHASTIC ERRORS; a -i 

- RANDOM DISTURBANCES AND SENSOR NOISE 

- TEST SIGNAL AMPLITUDE AND SHAPING 

- IDENTIFICATION TIME INTERVAL 
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STEADY-STATE IDENTIFIABILITY ANALYSIS (YARED) 

• Expected Likelihood Function 


I («) = E. 


K<">! 


= ^ [log det S+ tr(S~^S*)] 


• Residual Covariances 

S = E„ 

S. = E 


M" 


Kalman filter predicted residual covariance matrix 
= Actual residual covariance matrix 


Note: S» and !*(«) can only be computed when the true plant parameters are 

KNOWN. 
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EXPECTED MLE IDENTIFICATION SOLUTION 


• Expected Parameter Estimate (Theoretical) 

a. - E. {aj - ArG j"«N 

• Iterative Algorithms 

Gradient; 

Newton-Rhapson: °»j+i = “*j “ 


WHERE 


(a) = 
v2[*(a) 




STEADY STATE IDENTIFICATION ACCURACY 

• SYSTEMATIC ERRORS (BIASES) 

- PARAMETER ERRORS 

“bias = “ “* = 0 MODEL MISMATCH 

- INFORMATION MEASURE (YARED) 

I(a» ; a*) - I*(a*) - I*(a*) > 0 


• STOCHASTIC ERRORS 


I [-2, ■<;.>]-' E.{[% ci>] [|!t(»>]yi-<i..] 


[v2l*( a.)] 


-1 


(N+1) 


WHEN NO MODEL MISMATCH 



STOCHASTIC ERROR WITH PROCESS NOISE 



SIMPLIFIED IDENTIFICATION ACCURACY ANALYSIS 
(FOR EACH MODE) 


e[uu^] = UI 


E[w^] = WI 


. J 

k 

cb 

^ 



s^+2Cti)S+o)^ 



E[vv^] = RI 




a = (lu , 2Qni, b- , . . . b , c- , , , . f c^} - {b_ or c„} 

1 m' 1' p L M 


a 2^ 

U} 

= B?! 

1 

, 2.2 
(U) ) 

SNR 

(N+1) T 

“2^ 

4 

1 

(2Cto)2 

SNR 

(N+1) T 


1 

SNR 

[ m! 

a 2 
cm 

1 

SNR 



m 


A r I 

p = max { 


m j £ 1 

l,m |b|^ |c I ^ 


Ic|"[b|" U 




(MEAS. NOISE) 
(PROC. NOISE) 


(N+1) T 


(N+1) T 


HdE: THIS ANALYSIS ASSUMES THAT <-T « 1, c « 1 AND APPLIES FOR EACH MODE 




NUMERICAL RESULTS 


• PROBLEM SIZE 

12 MODES (=22 - 8R.B. - 2 SHALL) 

2 INPUTS s GIMBAL ANGLES) 

2 OUTPUTS (x> Y ATTITUDES) 

60 PARAMETERS (=12 MODES x 5 PARAMETERS/HODE) 




TEST SIGNAL. NOISE STATISTICS 
SAMPLE TIME: T = 0.1 SEC 

TEST SIGNAL! U = £|000(IN-LB)^ 

PROCESS NOISE; W = W(IN-LB)2 

MEASUREMENT NOISE: R = N x lO'^^ 


• WORST-CASE RELATIVE ERRORS AT TIME (MODE 9): 


PARAMETER 

1 SFC 

lA SEC 

(L39 HRS 


0.0265 

0.007 

0.0007 

2;oi 

3.75 

1.0 

0.1 

‘’1 

13.6 

3.6 

0.36 


2 . 7 

0.7 

0.07 

=2 

11.7 

3.1 

0.31 


SUMMARY 


• CONTROL PROBLEM 

- MODERN LQ CONTROL DESIGN WITH ROBUSTNESS RECOVERY 
PRODUCES ROBUST CONTROLLERS FOR LSS 


- ACCURATE ID ALLOWS A FIVE-FOLD INCREASE IN LOOP BW 


• IDENTIFICATION PROBLEM 

- STRUCTURAL MODES MAY BE IDENTIFIED ONE AT A TIME FOR 
SMALL DAMPING 


- LSS ID W/0 KF 

-GREATLY REDUCES PARAMETER BIASES 

-GIVES ONLY MODEST INCREASE IN STOCHASTIC ERRORS 

- RELATIVE ERRORS IN PARAMETERS AFTER ID ARE SMALLER FOR 
FREQUENCY THAN FOR DAMPING OR MODE SHAPES 

• OPEN ISSUE; HOW ACCURATE MUST ID BE FOR ROBUST CONTROL DESIGN? 
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THE DYNAMICS AND CONTROL OF 


LARGE FLEXIBLE SPACE STRUCTURES 


Peter M, Bainum, V. K. Kumar, R. Krishna, 
A. S. S. R. Reddy, and C. M. Diarra 
Howard University 
Washington, D, C. 


INTRODUCTION 


Large, flexible orbiting systems have been proposed for possible use in com^ 
munications, electronic orbital-based mail systems, and solar energy collection. The 
size and low weight-to-area ratio of such systems indicate that system flexibility is 
now the main consideration in the dynamics and control problem. For such large, 
flexible systems, both orientation and surface shape control will often be required. 

Figure 1 illustrates a conceptual development plan of a system software capa- 
bility for use in the analysis of the dynamics and control of large space structures 
technology (LSST) systems. This concept can be subdivided into four different 
stages: (1) system dynamics; (2) structural dynamics; (3) application of control 

algorithms; and (4) simulation of environmental disturbances. Modeling the system 
dynamics of such systems in orbit is the most fundamental component. 


SOLAR RADIATION PRESSURE EFFECTS 

The equations for determining the effects of solar radiation pressure on a 
flexible beam are summarized below. 

Forces : 


= -V / 




(t • n) ds (absorbing surface) 


" ^ ^ 2 

Fy = -2h^ / n(x • n) ds (reflecting surface) 


+ e(F^ - F^) 


(surface with reflectivity e) 


Moments : 


/X / 


(absorbing surface) 
s 


My = 2h^ J' n X r(t * n)^ ds (reflecting surface) 


M^ = M^ + e(My - M^) 


(surface with reflectivity e) 
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where 


T =s unit vector in the direction of solar radiation 
/\ 

n = unit vector normal to the surface ds 

h = 4.64 X 10“^ N/m^ 
o 

Results for a flexible beam; 




M„ = 


F,, = 


Y 




-hpa^bQ^Cz^ + Zq)!1 - 2/fljj |a(cos - cos + sin + sin 

/ 


(aoz’ - ^ r (a^z' - b^)^ 

-2h J z' dx i - 2h J ^ dx k 


(1 + z') 
2 


= 2h 


/• U„z'- V- / 

% (1 + z’)^ ^ 


where 

a = sin 6 b^ = cos 0 
o o 

z(x) = flexural deflection 



= nth modal frequency 
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Figure 2 shows the variation of the resultant horizontal and normal force com- 
ponents of a beam with a completely absorbing surface as the solar incidence angle 6 
is varied from 0^ to 90^. Here, 6 represents the angle between the normal to the 
undeflected beam and T. The horizontal and normal force components are measured 
relative to the beam^s undeflected axes. As expected, for small tip deflections of 
the beam, the resultant horizontal absorbing force component becomes zero for inci- 
dence angles of 0^ and 90^, while the normal component has a maximum amplitude at 
zero incidence angle. In figures 2 and 3 the individual effect of each mode is 
illustrated, with the assumed beam tip deflection as indicated in the figures. 

The magnitude of the resultant moments as the solar incidence angle is varied is 
shown in figure 3 for the assumed tip deflection of 0.01£. Large moments can result 
for larger deflections, whereas these moments would be zero for a rigid beam. For 
small pitch angle displacements, the moment due to solar radiation pressure may 
become greater than the moment due to the gravity-gradient forces, as shown in fig- 
ure 4. It is seen that at geosynchronous altitudes, the moment due to solar radia- 
tion may become predominant even for deflections on the order of 0.01£. With the aid 
of moment diagrams such as those in figure 3, it is possible to determine the dis- 
turbance torques due to solar radiation pressure once the number of modes and the 
associated modal deflections are specified for a model. 


MODELING ERRORS - ORBITAL AND GRAVITY- GRADIENT EFFECTS 

One of the principal sources of (disturbance) torques acting on an orbiting 
space structure is the orbital (gyroscopic) and gravity-gradient effects. Such 
effects associated with the orbital (angular) motion do not need to be considered 
when developing a system model for an Earth-based large flexible system. Many 
investigators, however, model the pitch, roll, and yaw modes (rigid body motion) of 
large, flexible orbiting systems as double integrator plants (two poles at the 
origin), and the subsequent control system design is based on these models. It is 
the purpose of this section to evaluate the effects of omitting the orbital and 
gravity-gradient effects when designing shape and orientation control laws for 
flexible systems in orbit. Models of flexible square plates and shallow spherical 
shells in orbit are selected as examples. 

The effects of designing control laws without the orbital and gravity-gradient 
torques included in the system models of square plates and shallow spherical shells 
in low Earth orbit (250 nautical miles) are illustrated in figures 5 to 10. A square 
plate was also considered in a geosynchronous orbit. The analysis was performed by 
first calculating the control law, which was in the form u = Fx, for the case where 
the orbital and gravity-gradient effects are not included in the model. The same 
control law is then applied to the model that includes these effects. 

For figure 4, the control law was selected such that the overall response time 
constant of the system is 2.22 hours (which may be reasonable for a large space 
structure). The shift in the closed loop poles of the plate model due to the pres- 
ence of the orbital and gravity-gradient effects is illustrated in this figure. It 
can be seen that some of the poles move to the right half plane, leading to insta- 
bility and thus emphasizing the importance of including the orbital and gravity- 
gradient effects in the model. The poles due to the rigid body modes are shifted 
considerably, but the flexible modes remain virtually unaffected. This can be 
attributed to the high frequency of the flexible modes. (Note that the orbital and 
gravity-gradient effects are of a relatively low frequency.) This result gives an 


94 



indication that by designing a more robust (faster response) control system, the 
shift of the actual closed loop poles would be relatively less pronounced. 

This phenomenon can be demonstrated by designing the gain matrix, F, such that 
the desired response time constant is reduced to 460 seconds. The shift in the 
closed loop poles for this case is shown in figure 6. In figure 7 the control forces 
are shown for the second (more robust) control law, where the closed loop response 
with orbital and gravity-gradient effects is degraded but does not become unstable. 

(It should be noted that time has been nondimensionalized with respect to the orbital 
frequency of the 250-nautical-mile low Earth orbit, in order to provide a basis for 
comparison.) The difference in the total control force impulse as applied to the two 
models (a) and (b) is minimal because of the robustness of the controller. The 
slowly varying orbital and gravity-gradient torques have a relatively greater impact 
on the less robust systems, and can even lead to possible instabilities, as was 
illustrated in figure 5. 

As expected, the orbital and gravity-gradient effects are less pronounced in the 
case of a structure in geosynchronous orbit than in the case of a structure in low 
Earth orbit. However, if the control systems are designed with response times com- 
parable to the orbital periods, under the influence of orbital and gravity-gradient 
effects the closed loop systems may become unstable. 

The shift in the closed loop poles of the spherical shell model due to the pres- 
ence of the orbital and gravity-gradient effects is shown in figure 8. One of the 
closed loop poles is moved to the right-hand side of the S-plane, causing instability. 
As compared to the case of the plate, the effect of the orbital and gravity-gradient 
torques on the shell is more pronounced, as the instability due to movement of the 
poles occurs at the relatively fast designed response time constant of 615 seconds 
(compared to 8000 seconds in the case of plate) . When the control is redesigned for 
a response time constant of 400 seconds, the shift in the poles is as shown in fig- 
ure 9. A general shift in the rigid body motion poles is observed, but the system 
remains stable. 

The control forces associated with both models (a) and (b) of figure 9 are com- 
pared in figure 10. A considerable increase in the control effort is observed when 
the model includes the effect of the orbital gyroscope and gravity-gradient torques. 
This may be explained by the fact that the mass distribution of the shell is more 
complex than that of the plate, resulting in relatively greater dynamic coupling when 
the gyroscopic and gravity-gradient effects are included in the shell model. 


THE DEVELOPMENT OF AN ALGORITHM TO EVALUATE COUPLING COEFFICIENTS 

FOR A LARGE FLEXIBLE ANTENNA 

The generic mode equations and the equations of rotational motion of a flexible 
orbiting body contain both coupling terms between the rigid and flexible modes and 
terms due to the coupling within the flexible modes that are assumed to be small and 
thus are usually neglected when a finite element analysis of the d 3 mamics of the 
system is undertaken. In this section a computational algorithm is developed which 
permits the evaluation of the coefficients in these coupling terms in the equations 
of motion as applied to a finite element model of a hoop/column antenna system 
(ref. 1). 
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Using a Newton-Euler approach, one can express the equations of motion of an 
elemental mass of the system, in the frame moving with the body, as: 

*|a^^ + r + 2w + r + a3Xr + 0 )x (Sx r)j^p dv = |f + e + L(q)/p|p dv (1) 

where 

p = mass per unit volume 
e = external forces per unit mass 

q = elastic transverse displacements of the element of volume 
f = force due to the gravity on the unit mass 

L = the linear operator which when applied to q yields the elastic forces 
acting on the element of volume considered 

r = position vector of element dv 

(A) = inertial angular velocity of the body frame 

== acceleration of the center of mass 
cm 


Equations of Rotational Motion 

The equations of rotational motion of the body are obtained by taking the 
moments of all the external, internal, and inertial forces acting on the body; l.e., 
from equation (1) : 


^r X + r + (2o3 x r) + (o) x r) + (o3 x (o) x r)^ p dv = ^ ^ ^ |l(q)/p + f + ^pdv 


( 2 ) 


One can obtain the following form for the equations of rotational motion: 


n=l n=l n=l 


( 3 ) 


where 


* y* jjo X (w X r^) - (r^ • io) (w x p dv 
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OO /• / .. 

X) “ J X q + 2r^ X (u X q) + X (w X q) + q X (u X 9^) 

V 

- (r^ • ui) Qai X q) - (q • (w X rQ>|p dv 


n=l 


oo 

= y* qp dv X (a 


n=l 


oo 

- f ) + 0 )^A A r X dv 

o' Z-rf n n y o 


n=l 


G^ = fr^ X Mr^ p dv 


°° r 

2 “ J (^o ^ ^ p 


n=l 


C = ^ r X e p dv 

V 


r = Tq + q 

M = matrix operator which when applied to r yields gravity-gradient forces 

a = acceleration of the center of mass 
cm 

f^ = force/mass due to gravity at the undeformed center of mass 

= modal shape vector for the nth mode 
60^ = frequency of the nth mode 

= time-dependent modal amplitude function 


Generic Mode Equations 

The generic mode equation is obtained by taking the modal components of all 
internal, external, and inertial forces acting on the body, i.e.. 




• • 


+ 2a)xr + 0)xr + a)x (ci)Xr) pdv 


>] 


= y* • jL(q)/p + f + Jp dv 


( 4 ) 
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The generic mode equation is obtained in the following form: 


\ + “X + + E = (sn + S 8nm + J An 

m=l \ m=l // 


where 


'^n ~ • w X ?o + • CO X (u X p dv 


OO 

= J' • 10 X q + • CO X q + • CO X (oj X q^ p dv 


m=l 


p dv 


( 5 ) 


S Smn = / 

m=l V 



V 


Here is the inertia coupling between rigid body modes and the nth structural 

mode and is the inertia coupling between the mth and nth structural modes. 


Cartesian Components of the Different Coupling Terms 

The expressions for R, G„, \p^, , g , and g in 

^ n mn n mn 

Cartesian components are presented in this section. The following vectors can be 
expressed in their Cartesian component form as: 


= ? i + 5 j 


+ C^k: 


(0 = CO i + 
X 


“y3 + 
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n=l 




V y ^ z 


X 


where i, j, and k are unit vectors along the body principal axes of inertia in 
the undeformed state, and ? , ? , and ? are the coordinates of a point in the 
undeformed state. ^ Y ^ 

With the use of the component forms of the vectors given above, one can expand 
the various vector expressions given in equations (3) and (5) to obtain 


R=[ja) +(J - J)coa)liH-rjd) +(J -1)03 0)^1 

Lxx y yzj Lyy ^ z^zxj-^ 

+ [j 0) 4* (J — J )(j0 o) ”1 k 

L z z y X X yj 


( 6 ) 


Xn) _ 


+ A - (h 

n \ yy zz / x \ : 



h(-)o) 

- H^^> 0 ) 


' zz 

/ X 

yx y 

ZX 



[(n) jj(n)V 

- 

+ 


xy 

yx / y 

\ ZX 

xz 

/ 

z 

1 

+ 0 ) 0 ) ( 

fiCn) ^ 

H(n)' 

y\ xz 

ZX } 

X z\ 

xy 

yx , 


+ - 0 ) 2 ) 

\ z yj\ yz zy j_ 


( 7 ) 




21 


( 8 ) 
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.(n) 


- *n[<»33 - «22>(“yf * “z”') ‘ * “zx ) 


t - i (h'"> + 

x\ yz zy 


+ 2M»_( 

31\ xy yx I 23\ 


+ M., ( + H 


+ - H 


23V“yy 


(n) 

zz 


)] 


) + 0) + (I) 

/ y\ zx xz / z\ xy yx / 


+ w to 
X y 


+ to to + to to (h^“^ + 

\^xy yx / y Z\ yz zy / z. x\ zx xz / 


( 9 ) 


x\ yy zz J y\ zz xx / z\ xx yy / 

tit = zkL + to + (0 

^mn ^%[jx\ yz zy ) y\ zx xz ) z\ xy yx ^ 

+ A fj) + J) + w 

m[”xV yz zy ) y\ zx xz ) z\ xy yx / 

+ to (0 + to to + to to 

+ “x y^xy yx / y z\ yz zy / z xV zx xz J 

- + lM- to^fL^”"^ + - to2(L<”“> + 

x\ yy zz / yV zz xx / z\ xx yy 


5„ - EE m 

a3 


where 


^mn " 




'6 


a3 J ^ot ^3 


a, 3 = X, y, z or 1, 2, 3 
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When a Is X in H 


(n) 

‘oi3 


when a Is 
for 3. 


z In H 


(n) 


a3 


or 

(mn) 

^3 ’ 

the 

a 

is y 

in 

or 

T (®n) 

a3 ’ 

a 


, the corresponding value of a in M „ Is 1. 


,(n) 
*a3 
Ls 3 


or 


(mn) 

^a3 ’ 


In M 


a3* 


a3 

a is 2 in and 

The same reasoning holds 


The expressions for and are obtained by the cyclic permutation 


of X, y, z In the expression for Q 


(n) 


In equation (7) , and the expressions for 


and are obtained by the cyclic permutation of x, y, z in the 

expression for G^ in equation (9) . 

For a discretized model, the expressions for the volume integrals are replaced 
by the following summations: 


“S’ - Z «a>l (fe”*)!”! 


" t Wii (♦S’)* 

i=l 


(a, 3 = x,y,z) 


( 12 ) 


(13) 


where 

k = total number of discrete masses 

i = index identifying a nodal point 

m. = mass concentrated at the ith node 
1 

£ = coordinates of m, in the undeformed state 

i 
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Figure !•- Development of system software for LSST dynamics analysis. 
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(b) Normal component . 


Figure 2.- Variation of solar force components with incidence angle. 

Totally absorbing surface - free-free beam (length Z = 100 m) . 
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Moment /width* M (N) 



Figure 3 


- Pitch moment due to solar radiation pressure 
surface). Effect of individual modes in the 


(completely absorbing 
system - free-free beam. 



Figure 4.- Moment due to gravity-gradient force as a function of 
pitch angle (100-m rigid beam) . 








(a) Designed closed loop poles. 
Model does not Include orbi- 
tal and gravity-gradient 
effects. 


Closed loop poles resulting 
from the control law of (a) 
when applied to a model which 
Includes orbital and gravity 


m 

Inag. axis 

gradient effects. 
9 

rv\\ 



O 27C 3TT5 T.o 


Real axis 


-9 

a 

-18 

1 ■ 

SI 

-27 

3) 

-33 

3 rigid + 3 flexible modes 

Figure 6.- Shift In closed loop poles due to orbital and gravity-gradient 
torques. Square plate In 250-n.-mi. orbit. Overall designed 
response time constant of the system = 460 sec. 




Forces 


fj, f*3 


Max. Force Astpllcudes 
f* - f - 14.60 N 

f 

2 - - 8.88 N 

- 22.50 N 
« 0.087 N 



(There is no appreciable difference In forces between 
model with G.G. and orbital effects and model without 
them.) 

ZZ f. (t) - 261.5 N-sec 
1 



(Mon-Dlmenslonallzed Time) 


Figure 7.- Time history of control forces. Square plate. 




>270 (a) Designed closed loop poles. 

Model doesn't Include orbital 
and gravity-gradient effects. 




(b) Closed loop poles resulting 
from the control law of (a) 
when applied to a model which 
Includes orbital and gravity- 
gradient effects. 
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1.0 


3 rigid + 6 flexible modes 



-180 


-270 


Figure 8.- Shift in closed loop poles due to orbital and gravity-gradient 
torques. Shallow spherical shell in 250-n.-ml. orbit. Overall 
designed response time constant of the system = 615.0 sec. 
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designed response time constant of the system = 400 sec. 
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PEAK FORCE AMPLITI]DES-C0NTR0L LAW DEVELOPED WITHODT 
ORBITAL AMD GRAVITY' GRADIENT EFFECTS IN THE MODEL AMD 
THEM APPLIED TO A MODEL WITH THEM (Shell in orbit). 


Without Orbital 
and G.G. Effects 

With 

6« 6* 

Orbital and 
Effects* 


'r 

516.40 N 

'i*= 

565.63 N 



73.30 N 

*2*= 

164.76 N 



239.50 N 

'3*= 

321.11 N 


*4= 

117.54 N 


117,07 N 


*5= 

f,: 

132.45 N 


153.84 N 


146.82 N 

431.61 N 


7020.8 N-sec 


a n 


19373 ,0 N-sec 


T (• a)t) 

c 

(Non-Dlmenslonalized Time) 


Figure 10.- Time history of control forces. Shallow spherical shell. 




CONTROL OF STRUCTURES IN SPACE 


L. Meirovitch and H- Barun 
Department of Engineering Science and Mechanics 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 


EQUATION OF MOTION FOR DISTRIBUTED SYSTEMS 


Differential equation: Lu(P^t) + M(P)3^u(P^t)/Jt^ = f(P>t) ^ Pe D 

Boundary conditions: Bju(P>t) = 0^ l=L2>...^p ^ Pes 

u(Pyt) = displacement at point P 

Ij Bj = differential operators (L is self-adjoint of order 2p) 

M = mass density 

f(P^t) = distributed control force 


EIGENVALUE PROBLEM 

Differential equation: L4*=AM+ 

Boundary conditions: Bj«#»= 0, i = L2^...,p 
Solution: eigenvalues A ^ =o^ , eigenfunctions (r=L2.,,.) 
Because L is self-adjoint^ eigenfunctions are orthogonal 
L is generally positive semldefinite -► Ap are all nonnegative 

Ap = 0 for rigid-body modes 

= natural frequencies 

Orthonormality conditions: M<f»s+r clD = Jpg . fgLfp dD = Apfp^ 
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MODAL EQUATIONS 


Expansion theorem: u(P>t) = Z 

r=i I I 

Up(t) = modal coordinates 

Modal equations: Up(t) +o^Up(t) = fp(t) > r=L2>... 

Modal controls: fp(t) +p(P)f (Pyt)dD ^ r=h2^... 

Coupled controls: fp(t) = fpiu^^ Ui^ U2-' ■' f=1^2^... 

Independent modal-space control (IMSC): fp(t) = fp(Up^Up)^ r=L2^... 


CONTROL IMPLEMENTATION 

Distributed actuators: f(Pyt) = Z 

Coupled controls: unable to design distributed controls 

IMSC: design modal controls firsts then use above formula 

no control spillover 

m . 

Discrete actuators: f(P>t) = Z F|(t) AIP-P.) 

j = l J J 

m m 

f.(t) = Z = 21 BriFj(t) > r=L2..... 

r j=i T'r j j j=i rj j 

n = number of controlled modes 


i-[Pl F2 




f = 0, 


fjl B - B 


ri 


Coupled controls: Design f(t) so as to ensure controllability 
IMSC: Design f (t) so^as to control a given number of modes 

Then^ FTt) = B^f(t) . B^ = pseudo-inverse of B 
To avoid pseudo-inverses^ let m = n^ or the number of actuators 
must equal the number of controlled modes. 


CONTROL IMPLEMENTATION (CONT^D) 


Distributed measurements: 

Measurements: u(P^t)^ u(P/t) for all P and at any t 

Then^ modal coordinates and velocities^ Up(t) and Up(t)^ are computed 

by using the modol filters 

Ur(t) = ( H(P) +r(P)U(P.t)dD. Ur(t) = ( M(P) r=L2. . . . 

Discrete measurements: 

Measurements: yj(t) = u(Pj^t)^ yj(t) = u(Pj^t)^ j = h2>...A 
k = number of sensors 

Standard approach: Use Luenberger observer to estimate state 

Discrete measurements treated as distributed: 

Use interpolation functions to compute estimate ti(P/t) of u(P^t) 

Then^ use modal filters to compute estimates Uj-(t) of Up(t) 

Divide structure into s segments (elements) 

A ^ T 

Approximate displacement: u(P^t) = E L (P)Vi(t) 

j=l 

Vj = measurements at the boundaries of j'th interval 

Lj = vector of interpolation functions (from the finite element method) 

Estimated modal coordinates: 

s s 

Ur(t) = S M(P)+r(P) H L|vi(t)dD= E IriVi(t). 
r r -j~j -rj-j 

r = ly2y...yn 

c s 

Similarly Up(t) = E IriVj(t) 
r j = j ^ f J ^ J 

irj " lj^M(P)+p(P)Lj(P)dD = const 

are computed off-1 ine^ in advance. 
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CONTROL IMPLEMENTATION (CQNT^D) 


Rearrange Iri such that 

Up(t) = i Cfjyj(t) . Uf-(t) = ^rJ 

Let 

u(t) = [uj U2‘" ^2 • • ■ Up]^ 

C “ C p j J r~ 1 ^ 2 ^ ■ • » j rij j~ 1 ^ 2 ^ ■ • i^k 

y(t) = [yi y2 ... y(t) = [vi 1/2 .•■ yj^ 

u(t) = Cy(t) y u(t) = Cy(t) 

M 

The way C is assembled depends on the nature of the interpolation 
functions (see example later). 


THE LANGLEY BEAM EXPERIMENT 


Free-free beam controlled by using 4 actuators and 9 sensors 

Use IMSC to control four modes^ two rigid-body and two elastic modes 

Actuator forces Fj(t) = (B"^)f-jff(t)^ j = L2^3>4 

Modal forces: 

1) For rigid body modes 

7 p = I Upl +1 Upl/Cp . r=h2 
Cp = weighting factor 
if 7 p < dpy then fp = 0 
if 7 p > dp and 

i) Up > 0> Up > 0> or Up > 0 > Up and lUpl< fp> then fp = -kp 

i i ) Up < 0 ^ Up < 0 ^ or Up < 0 < Up and I Up)< c^j then f p = kp 
dp = magnitude of the deadband region 

£p = threshold velocity^ kp = modal control force 



THE LANGLEY BEAM EXPERIMENT (CONT'D) 


2) for elastic modes 

fj-(t) = -kj. > |Ufl » dp and Up > 0 

0 y -dp < Up < dp 

kp ^ / Up I ^ dp and Up < 0 

For simulation purposes^ the response is available in closed form both 
for rigid-body and elastic modes. 

Actuator locations immaterial when IMSC is used. 


F(t) = B"^f(t) Fj(t) = z: (B‘^irfr(t) j=L2,3.4 
~ ~ J r=l ' 


Actual control forces are a combination of on-off functions of the type 


-Jr 


fr 

Ur 

Jr 7r 

•Ur 


Actuators available at LARC have four components/ which can be assigned 
to each modal control force. As a result/ each actuator command becomes 
a linear combination of 4 modal on-off control forces. 



THE LANGLEY BEAN EXPERIMENT (CONT'D) 


Sensors measure displacements alone. Ideal ly> for bending they should 
measure displacements and slopes (velocities and slopes of velocities too) 

Divide the beam into four equal segments (elements) and measure displacements 
at the ends and middle point of the elements so that the nine sensors are spaced 
at equal intervals. 

As interpolation functions use 

Lj = [i(2^D l-3t+ 0< t< \ 

where ^ is a local coordinate related to the global coordinate x by 

i = j - in which j is the segment number 


The C matrix is assembled fromJ^pj tensor as 

^rp " Sl-rJi ^2j+l-2>p ^ ! U2,...jn , n=A 

p = \j2j k = 9 


£ = index denoting the interpolation function 
j = index denoting the element number 

Because velocity measurements are not available estimate velocities by using the 
relation 

Ur(jT) - Ur(JT-T) 

Up(jT) = 

T = sampling time = 1/33 sec. 

Or^ one could use a modal Luenberger type observer. Because the controls 
are nonlinear^ the convergence of the observer can only be determined by trial 
and error. 
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THE LANGLEY BEAM EXPERIMENT (CONT*D) 


Parameters associated with the beam: 

L = 12 ft. cross-section = 6 x 3/16 in 

6061 aluminum: p= 0.1 Ib/in^ E = 1 x 10^ Ib/in^ 


The free-free> 
transcendental 

Oj = 0 
^2 = 0 

0 J 3 = 11.47979 
= 31.64450 
W 5 = 62.03586 
<Jg = 102.5484 
«7 = 153.1897 


uniform beam admits a closed-form solution, 
equation was solved numerically to yield the 


rad/s 

rad/s 

rad/s 

rad/s 

rad/s 


The 

eigenvalues 


Simulation of the Beam motion : 

The first 7 modes are Included in the simulation: 4 controlled + 3 
residual modes 

Control Gain Parameters: 

dj = d 2 = 0.002 . d 3 = di( = 0.0005 

kj = l <2 = 0.3 . k 3 = 0.12. ki, = 0.03 

= £-2 = O'Ol 

= lUpl + lu^l /lO . r=1.2 

Sampling time = 1/33 sec. 

Viscous damping was added to each flexible mode 
Damping factor = 0.002 . r=3.4.5.6.7 

Disturbance of the beam was taken in the form of a unit Impulse of magnitude 1/12 lb 
applied at Xg = 0.67L. 

The displacement of the beam cannot exceed 1 in at any time because of the location of 
the sensors and actuators in the experimental setup. 
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THE LANGLEY BEAM EXPERIMENT (CONT*D) 


Results 


• The main contribution to the response is from the rigid-body modes. 

•The second elastic mode shows noticeable participation. This participation will 
eventually disappear due to internal damping. 

• Control of the second elastic mode can be enhanced by sensing velocities^ or 
estimating velocitites via a Luenberger observer. 


• Observation spillover (which may arise from the need of more sensors) was found 
to be negligible. So was the control spillover into the residual modes. 
Simulations of the beam with and without the residual modes indicated that 
spillover effects are infinitesimal. The reasons for this are: 

1) IMSC is used 

2) Modal filters are used 

3) Residual modes have very high frequencies 


Conclusion 

The IMSC method in conjunction with on-off modal controls does a good Job 
in controlling the motion of the beam^ where the motion consists primarily of the 
rigid-body modes. 
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ROBUST PRECISION POINTING CONTROL OF 
LARGE SPACE PLATFORM PAYLOADS 


S. M. Joshi 

Old Dominion University Research Foundation 
Norfolk, Virginia 



LSS CONTROLLER DESIGN STRATEGY 


INHERENT DAMPING PLAYS A VERY IMPORTANT ROLE IN STABILITY, 
THEREFORE IT IS HIGHLY DESIRABLE TO ENHANCE INHERENT DAMPING 
USING A SECONDARY CONTROLLER. 

DESIGN PRIMARY CONTROLLER FOR CONTROLLING RIGID-BODY MODES 
AND POSSIBLY SOME SELECTED STRUCTURAL MODES. 


SECONDARY CONTROLLER - VELOCITY FEEDBACK 

flexible part of DYNAMICS: 

D = ? O ^ A. = , ' ' ' ) 

SENSOR output: 

>{• = f9 

CONTROL LAW; 

/■= = -^^9 

WITH PERFECT (LINEAR^ INSTANTANEOUS) ACTUATORS/SENSORS : 

- STABLE IF ^ O 

- ASYMPTOTICALLY STABLE (AS) IF 

K), AND^1/1^^^ CONTROLLABLE 

WHAT IS THE EFFECT OF 


NONLINEARITIES 

actuator/sensor dynamics 



EFFECT OF SENSOR/ACTUATOR NONLINEARITIES: 



ACTUATOR nonlinearities : LET ^0^ = O [ ^ ACTUATOR j 

a) ORIGIN STABLE IF 2 O 

b) ORIGIN ASIL* IF K"r 

AND CONTROLLABLE 

(for DIAGONAL %i^.) > o J 

NONLINEARITIES IN ACTUATORS AND SENSORS : 

LET Ky. BE DIAGONAL 

a) ORIGIN STABLE IF C" ^(<T) > 0 FOR EACH ACTUATOR AND SENSOR 

b) ORIGIN ASIL* IF (T " " 

AND (-A.,0^) IS CONTROLLABLE 

CONSIDER actuator/sensor DYNAMICS GIVEN BY 

~ -^a Bck^a 

WHERE yU >0, SMALLJ STRICTLY HURWITZ 

• IF THE CLOSED-LOOP SYSTEM WITH PERFECT ACTUATORS/SENSORS 

IS AS^ THEN THE TRAJECTORY WITH FINITE-BANDWIDTH ACTUATORS/ 
SENSORS IS OC/0 CLOSE TO TRAJECTORY WITH PERFECT ACTUATORS/ 
SENSORS (ref. 1). 

HOW TO DETERMINE /U THAT WILL GUARANTEE STABILITY? 

* ASYMPTOTICALLY STABLE IN THE LARGE 
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• CONSIDER THE SINGLE-INPUT SINGLE-OUTPUT CASE (SISO) 



ASSUME SIMPLE STRUCTURAL MODE FREQUENCIES^ NO DAMPING^ AND 
NO POLE-ZERO CANCELLATIONS (i.E. CONTROLLABILITY) 

• THEOREM the CLOSED-LOOP SYSTEM IS AS FOR SUFFICIENTLY 
SMALL >0 IF 

AT C < = i, ‘ 

WHERE = /(SaO^) 

• ALWAYS STABLE FOR 1ST ORDER GqCs) 

• when actuator/sensor have Hci real poles 

AT 5=-(5a and no zeros, the system is as for 

SMALL Ky.>OlF 

oi, > 

fan Cnyzna') 

GIVES some insight INTO ACTUATOR/sENSOR BANDWIDTH 
REQUIREMENTS. ADDITIONAL INVESTIGATION NEEDED FOR 
OBTAINING MORE USEFUL RESULTS. 

(MIN. REQD.) 
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CONTROL OF PPS/LSS 


• ATTITUDE AND VIBRATION CONTROL OF LSS 

• POINTING CONTROL OF EACH PPS 

• SINCE THE MASS OF EACH PPS MAY BE OF THE SAME 
ORDER AS THAT OF LSS^ INSTABILITY IS POSSIBLE 
IF CONTROL SYSTEMS ARE DESIGNED INDEPENDENTLY. 


SECONDARY CONTROLLER USING ANNULAR MOMENTUM CONTROL 
DEVICES (AMCD'S) 


ASSUMPTIONS 

• AMCD RIMS ARE RIGID 

• RIM DIA = 2 M (small COMPARED TO LSS) 

• ACTUATORS AND SENSORS PERFECT 

• ELECTROMAGNETIC ACTUATORS AND POSITION SENSORS 
ARE ALMOST PERFECTLY LINEAR IN THE OPERATING 
RANGE. BANDWIDTH IS SEVERAL HUNDRED Hz. 

• ACTUATORS/SENSORS COLLOCATED 

• INHERENT DESIGN CHARACTERISTIC OF AMCD 

• CONTROL law: 

WHERE S IS THE ACTUATOR CENTERING ERROR VECTOR 

• STABLE IF Kp >0 , Kj- ^ ® 


ASYMPTOTICALLY STABLE (AS) IF 


b) LSS STRUCTURAL MODEL STAB I LI ZABLE 

c) NO UNDAMPED LSS MODES AT ^j-^l = SPJN FREQ, 

V 

d) ^ ^ 

I 

PRIMARY ATTITUDE CONTROLLER (PAC) 

• USING TORQUE ACTUATORS 

• COLLOCATED ACTUATORS/SENSORS 

• NONCOLLOCATED ACTUATORS/SENSORS 

• USING AMCD'S 


PAC USING TORQUE ACTUATORS 

• SEVERAL TORQUE ACTUATORS AT VARIOUS POINTS OF LSS 

• COLLOCATED ACTUATORS/SENSORS : 

• EQUATIONS OF MOTION! 

Qx, = rr 

• COLLOCATED ACTUATORS/SENSORS 

MEASUREMENTS CONSIST OF: 

• ATTITUDE VECTOR: ci ~ P 

• ATTITUDE RATE VECTOR: o( Z 

CONTROL LAW: T= - f <x' <3,. *9 
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• STABLE IF Gp , Gf 2 O 

• AS IF Gt.?0, (C^ ^ P”^ CONTROLLABLE 

• CONTROL LAW MINIMIZES 

<50 

; rn 

0 



o 


o 

Qt 

o 


s 

o 


/-X5 

^5 


At 


• EFFECT OF ACTUATOR/sENSOR DYNAMICS 

SISO CASE^ NO REPEATED FREQS. ^ CONTROLLABLE 

^et , ^a(o)=0 

Gp - f^Jp , 

• AS FOR SMALL k ~>0 IF 

- <Pix ^ < iSo- 4 . 

WHERE ^ FUNCTION 

OF (jOj^ AND POSITION RATE GAINS 


• ADDITIONAL INVESTIGATION IS NEEDED 


• IF TORQUE ACTUATORS AND ATTITUDE/RATE SENSORS ARE 
NOT COLLOCATED 

• STABILITY NOT GUARANTEED 

• MUST USE LQG-BASED APPROACHES INVESTIGATED EARLIER 

• KNOWLEDGE OF FREQUENCIES AND MODE-SHAPES REQUIRED 

A TRUNCATION: "rESIDUAL" MODES IGNORED 
IN THE DESIGN PROCESS 
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A MODEL ERROR SENSITIVITY SUPPRESSION (REF. 2 ), 
*SPILLOVER* IS INCLUDED IN THE QUADRATIC 
PERFORMANCE FUNCTION: 


= Ac 

= AyXr 





• IT CAN BE SHOWN THAT tHE SAME AMCD'S CAN BE USED FOR 
PRIMARY CONTROL 

• ATTITUDE AND RATE SENSORS LOCATED ON LSS AT MIDPOINTS 
OF AMCD ACTUATORS 

• APPROXIMATES TORQUE ACTUATORS AND COLLOCATED SENSORS- 
PROVEN TO BE STABLE 

• CLOSED-LOOP RIGID-BODY BANDWIDTH DEPENDS ON TOTAL 
MOMENTUM^ ALLOWABLE GAPSj ETC. 


PRELIMINARY HATH MODEL OF LSS/MPPS 


• ASSUME LUMPED POINT-MASSES AT POINTS OF ATTACHMENT OF 

PPS'S (for LSS model) 

• ALL PPS RIGID 

• EACH PPS CONSISTS OF GIMBALS AND TORQUERS (ONLY X-AXIS 
GIMBAL ASSUMED IN PRELIMINARY STUDY). 
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• MATH MODEL GIVEN BY: 





Ax + 3x t Cx ~ V f 


^ n<yx I 



= diAj col, • • '* ^ 

f = ‘ .rj, 

^ - (-Txi ^Tyi)^ 


OBSERVATIONS 

• INERTIAL ATTITUDE AND RATE SENSORS ON EACH PPS PAYLOAD 

• » INERTIAL ATTITUDE AND RATE SENSORS ON LSS 

(collocated with LSS torque actuators) 

• SENSORS FOR MEASURING RELATIVE ANGLE BETWEEN EACH 
GIMBAL AND LSS 


CONTROL OF MULTIPLE PRECISION - POINTED STRUCTURES (MPPS) 

MOUNTED ON LSS 

• MPPS USED FOR COMMUNICATIONS^ ASTRONOMY^ EARTH RESOURCES, 
AND WEATHER PAYLOADS 
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advantages: 


• LESS ORBITAL SLOT SPACE AT GEO 

• SAVINGS IN POWER^ CYROGENICSy GROUND DATA LINKS^ AND 
SMALLER GROUND TERMINALS 


CONTROL LAWS 

CONTROL LAW I: DECENTRALIZED CONTROL OF LSS/PPS: 

• USE LSS ATTITUDE AND RATE SIGNALS AND DESIGN LSS 
CONTROL LAW (COLLOCATED ACTUATORS/SENSORS) 

• USE PPS ATTITUDE AND RATE SIGNALS AND DESIGN 

CONTROL LAW (TO GENERATE DESIRED GIMBAL TORQUER TORQUES) 

• THE RESULTING CLOSED-LOOP SYSTEM CAN BE UNSTABLE 

CONTROL LAW II: ROBUST COMPOSITE CONTROL 

• USE INERTIAL PPS AND LSS SENSORS^ AND ALSO GIMBAL- 
ANGLE SENSORS AND COMBINE THE SIGNALS TO OBTAIN 

y - r’x 

• CONTROL law: 

f = - Kpy 

• CLOSED-LOOP SYSTEM Lyapunov-stable // 

• ASYMPTOTICALLY STABLE (AS) IF f^p,Ky 

LSS STRUCTURAL MODEL STABILIZABLE^ REGARDLESS OF 
NUMBER OF MODES AND NUMBER OF PPS. THEREFORE^ 

CONTROLLER IS ROBUST. 



NUMERICAL RESULTS 


100' X 100' X 0.1" COMPLETELY FREE ALUMINUM PLATE 
TWO PPS EACH WITH MASS “ LSS MASS 

THREE LSS TORQUE ACTUATORS WITH COLLOCATED ATTITUDE AND 
RATE SENSORS 

CONTROL LAW I: DECENTRALIZED CONTROL 

• LSS BW “0.05 raa/sec 

• PPS‘ BW INCREASED GRADUALLY - A STRUCTURAL MODE WAS 

DRIVEN UNSTABLE FOR PPS BW > 0.1 rad/sec 

CONTROL LAW II: COMPOSITE CONTROL 

• LSS BW “ 0.05 rad/sec 

• PPS BW OF 1 rad/sec 0.7) 

WAS EASILY OBTAINED WITHOUT SIGNIFICANT 
EFFECT ON CLOSED-LOOP LSS STRUCTURAL MODES 


CONCLUDING REMARKS 

TWO-LEVEL LSS CONTROL IS STABLE AND ROBUST AND 
OFFERS PROMISE 

FURTHER INVESTIGATION NEEDED ON EFFECTS OF 

actuator/sensor bandwidth 


COMPOSITE LSS/MPPS controller is stable and ROBUST 



PLANS FOR CONTINUED RESEARCH 


• COMPLETE INVESTIGATION OF LSS/HPPS COMPOSITE CONTROL^ 
INCLUDING PERFORMANCE EVALUATION 

• INVESTIGATE ANNULAR SUSPENSION AND POINTING SYSTEM (ASPS) 
FOR PPS CONTROL ACTUATION 


• START INVESTIGATION OF HOOP-COLUMN ANTENNA CONTROL 
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Magnetic 
actuators Z 



AMCD/LSS CONFIGURATION 
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IMAGINARY 




NONLINEAR EFFECT ON MODAL DATA ANALYSIS METHOD 


Lucas G. Horta 

NASA Langley Research Center 
Hampton , V ir ginia 



NONLINEAR EFFECT ON MODAL DATA ANALYSIS METHODS 


OBJECTIVES: 


1 - DETERMINE HOW THE PRESENCE OF NONLINEARITIES IN STRUCTURAL 

TEST DATA CAN BE DETECTED WHEN USING MODERN LINEAR MODAL DATA 
ANALYSIS METHODS. 

2 - EVALUATE THE EXTENT TO WHICH LINEAR ALGORITHMS CAN PROVIDE 

USEFUL INFORMATION ON NONLINEAR SYSTEMS. 


APPROACH: 


1 - GENERATE SIMULATED TEST DATA BY A NONLINEAR ANALYTICAL MODEL. 

2 - USE LINEAR METHODS (IBRAHIM TIME DOMAIN ALGORITHM (ITD) AND 

FREQUENCY-DOMAIN TRANSFER FUNCTION TECHNIQUES) TO ANALYZE SETS 
OF THIS DATA WITH CONTROLLED PARAMETRIC VARIATION. 


SINGLE DEGREE OF FREEDOM MODEL 



S'Y&TEM Pi^^ePfeMTiM £C|^uMrON •foft Pt=>p\ftC6 



(aeMePAl Form 

e« c, 9® + Ci6 + c^e 

FRCQueMCy S 0.2.89 

LinjEaQ ModAl DamPiN 6 a 0-01 
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SINGLE DEGREE OF FREEDOM RESPONSE 



TIME CSEO 

C'/Ca"- 0 



TIME (5EC) 


C>/C2= 10 


THE IBRAHIM-TIME-DOMAIN (ITD) MODAL DATA ANALYSIS METHOD 


0 AN OFF-LINE LARGE-SCALE DATA ANALYSIS METHOD DEVELOPED FOR STRUCTURAL 
DYNAMICS TESTS. 


0 OPERATES ON FREE-DECAY RESPONSES SOLVING MANY DATA CHANNELS SIMULTANEOUSLY. 


0 HAS PROVEN MORE SUCCESSFUL THAN OTHER LABORATORY METHODS FOR HANDLING NOISE, 
LARGE SYSTEMS, AND CLOSELY SPACED MODES. 
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SINGLE DEGREE OF FREEDOM RESULTS 



C1/C2 = 5 


ITD 


12 

10 


© 


IDENTIFIED 

DAMPING. 

% CRITICAL 



o 

0 

o 


10 15 

© 


INIT. AMPLITUDE (DEG.) 

C1/C2 = 5 


SINGLE DEGREE OF FREEDOM RESULTS COMPARISON 


' 

Fl (Hi) 

F*2 { hO 

Fi (hz) 

Ana lyTicAl 

0.104Z 

0.^127 

I.5Z 

ITD 

0.^030 

O.8S70 

\-5<ol 

^difference 

o.i<^ 

1.72 

2.7^ 
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TWO DEGREE OF FREEDOM ANALYTICAL MODEL 


Syi>TEM 


^ \ 




1 

n 



I] 





\ 'S. ^ 'V 'S. "V 



L)^ 4 E^R Mode 
- . 100 Wtl 


»0M of MoTiOKi 



\ \ V V X 



^ 

NJon liNicAR Mode 
Pz“ 


TWO DEGREE OF FREEDOM RESULTS 



Assort ED DOF = Co y\ SSUmED DOF - lO 


tP3M= EMER-CaV pROviPED To bECOMD M\OD£ 


139 



SUMMARY 


1 - THE ITD MODAL DATA ANALYSIS METHOD SUCCESSFULLY IDENTIFIED THE FREQUENCY 

COMPONENTS (BUT NOT DAMPING) OF THE TRUE SOLUTION OF A NONLINEAR SYSTEM 
FROM SIMULATED TEST DATA. 

2 - THE APPROXIMATE MODAL PARAMETERS OF THE LINEAR SYSTEM CAN BE IDENTIFIED BY THE 

ITD METHOD WHEN LOW LEVELS OF NONLINEARITIES ARE PRESENT. 

5 - THE ITD METHOD SUCCESSFULLY IDENTIFIED THE MOD-AL PARAMETERS OF A LINEAR 
MODE IN THE PRESENCE OF A NONLINEAR RESPONSE. 


FUTURE RESEARCH THRUST 


1 - EXPERIMENTAL EVALUATION OF ITD METHOD ON A TW0-DE6REE-0F-FREED0M NONLINEAR 

LABORATORY MODEL. 

2 - APPLICATION OF FREQUENCY-DOMAIN TRANSFER-FUNCTION TECHNIQUES TO THE SAME 

ANALYTICAL AND EXPERIMENTAL DATA. 


STRUCTURAL DESIGN FOR DYNAMIC RESPONSE REDUCTION 


Brantley R. Hanks 
NASA Langley Research Center 
Hampton, Virginia 



OBJECTIVE: STUDY STIFFNESS AUGMENTATION BY MATHEMATICAL DESIGN 

APPROACH: APPLY LINEAR REGULATOR THEORY WITH PROPORTIONAL FEEDBACK 

JUSTIFICATION: STIFFNESS IS READILY AVAILABLE TO DESIGNER AS PREDICTABLE 
PASSIVE CONTROL 

TIME-INVARIANT LINEAR REGULATOR— GENERAL 


SYSTEM: 

X = Ax + Bu + Dw 

CONTROLLED VARIABLES: 
y = Cx 
OBJECTIVE: 

Min J where J = 
u 

OPTIMAL CONTROL (ASSUMING w IS RANDOM!): 
u = -R"^ b"*^ Px 
WHERE P IS SOLUTION TO 

P = -PA-A^P + PBR"^ B^P - c"^QC PCtf) = 

IF GET STEADY-STATE P (AND U) FROM 

0 = -PA - A^P + PBR"^ B^P - C^ QC 

POSITIVE DEFINITE P EXISTS IF 

• A IS DETECTABLE IN C, STABILIZABLE IN B 

• RESPONSE WEIGHTING MATRIX, Q, IS POSITIVE SEMIDEFINITE 

• CONTROL WEIGHTING MATRIX, R, IS POSITIVE DEFINITE 


T Sj Xj 


+ ( [y*'Qy + Ru j dt 



LINEAR REGULATOR ADAPTED TO STRUCTURES 


SYSTEM; 


f 

X 

W 


-M -M 


X 


X 


> + 




u 


M 


w w ~ N( 


OBJECTIVE FUNCTION: 

I 

J = 

ASSUME: 


•"x 

1 

K 

0 

f.1 + (uT)[e]W' 

L 

^ 0 

M __ 



dt 


• RANDOM INITIAL CONDITIONS 

• COMPLETE STATE FEEDBACK WITH NO ROTATIONAL COUPLING 

• t^— ► » (TIME-INVARIANT STRUCTURAL CHANGE) 


CONTROL: 



WHERE P^ ^ AND 


0 0 
-m'^br’^bT m~^ -m'^br'^b"’’ m“^ 

P 22 are SOLUTIONS TO 



pT A + aT pT _ p M ^ 
21 21 21 21 21 ^ 


BR ^b"*^(m ^)^P23^ = 0 


AND 


^ 22^22 ^ 22^22 


P 22 M ^BR ^ 


bT(m 


P + ( P 
22 ^ 21 


+ P 


21 


C 2 Q 2 C 2 ) 


IN THESE EQUATIONS A^j^ = M ^ K AND = M ^ G 

NOTE THAT (1) IS NOT SYMMETRIC; ALSO THAT (1) IS INDEPENDENT OF (2). 



Q WEIGHTING MATRIX CONSIDERATIONS 



• If rate and displacement considered independently and Q chosen so as not 
to couple X and x 



• For design , selection of C is governed by desired minimum response points. 
Hence, C and Q may be assigned similar functions. 

• Diagonal C and Q minimizes weighted square response at selected coordinates. 

• Choice of Q = K and = M minimizes sum of strain and kinetic energy at 

n zz 

locations determined and (optionally) weighted by C. 
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REGULATOR FOR STRUCTURES-MODAL COORDINATES 


TKANSFOEMATION x = (j>q WHERE q = q e ' 

K 

WHERE (J. IS NORMALIZED (fi'*M(|) = I 

AND IS ASSUMED PROPORTIONAL TO (i.E., = -25^^01^ OR G(|) 

OBJECTIVE FUNCTION BECOMES 


J 


NOTE THAT 



HENCE, 



WEIGHTING MATRIX Q = 



0 




RICCATI EQUATIONS BECOME 




-c\b^C^ - 0 ( 3 ) 

AND 

'’32^50.] + 3- ^'■'"^'’22 - ('’21 ■ " W 

WHERE P , B, R, and C aRE MODAL EQUIVALENTS OF P, B, R, AND C. 

BY CHOOSING P, B, R, AND C DIAGONAL, WE DECOUPLE THE SOLUTION AND GET PURE "MODAL 
CONTROL." 
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CANTILEVER BEAM MODEL 





J xj. I 


|x^ 


ASSUMED: 

• CONSISTENT MASS FINITE ELEMENTS 

• UNIFORM INITIAL STIFFNESS & MASS DISTRIBUTION 

• FIRST NATURAL FREQUENCY = .OA7 Hz (.297 rad/sec) 


PHYSICAL IMPLEMENTATION OF STIFFNESS CONTROL 


COUVUNO 


r 





TENSIONED CABLE 




•COMPRESSION MEMBER 


A COUPliNG- 






1 : “ 

J 

/ 




PRACTICAL OPTIMUM ? 
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CONTROL WEIGHTING EFFECTS ON DESIGN 


UNDAMPED NATURAL FREQUENCIES 



INITIAL 


FINAL FREQUENCY, 

SEC 

MODE 

FREQ,,™ 

R = loi 

R = I 

R = .1 I 

R = .01 I 

1 

.297 

.359 

.557 

.972 

1.725 

2 

1.867 

1.880 

1.989 

2.619 

4.538 

3 

5.262 

5.267 

5.309 

5.684 

7.711 

4 

10.382 

10.384 

10.406 

10.615 

12.233 


DAMPING RATIOS 


MODE 

INITIAL 
DAMPING 
% C/C^ 


FINAL DAMPING, Z 

'CR 

R = 101 

R = I 

R = .1 I 

R = .01 I 

1 

2 

59.3 

108 

176 

298 

2 

2 

12.0 

35.0 

78.1 

131 

3 

2 

4.7 

13.4 

38.6 

82.9 

4 

2 

2.9 

7.1 

21.0 

55.4 


•note: solutions obtained separately for stiffness and damping 

COMPARED EXACTLY TO FULL ORDER CONTROLLER SOLUTION 


147 


STIFFNESS MATRIX COMPARISON (ASSUMED CONSTANT MASS) 


ORIGINAL K 


X| 02 X 2 02 

125 -1250 -125 -1250 

16667 1250 


X3 03 


250 


0 
I 

8333 i_ _0 
0 -125 

33333 1250 
250 


Xz 

0 

0 

0 

0 


0 / 


0 

0 

0 

0 


0 
0 

1 

-1250 ' 

I 

8333 I 
0 -125 -1250 

33333 1250 8333 

250 0 

33333 


FINAL K FOR R = .1 I 


L 


-1261 

-127.8 

-ml -1.7^1 
1 

-3.83 

-.166 

-.71" 


16719 

125A 

8278 1 6. ill 

1 

12.85 

.84 

2.94 



260.3 

- 5.22 -127 

-1239 

“I 

1 -1.58 

-3.84 

X 2 



33528 1236 

82A9 

1 2.45 
1 

4.58 

02 



2A6.3 

-2.06 

-124.6 

-1239 

^3 

Ui - , 


33540 

1234 

8252 

% 





268 

-3.37 

X, 






33543 

9., 
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R WEIGHTING EFFECT ON STIFFNESS MATRIX 
(FIRST ROW ONLY SHOWN) 



RELATED SPONSORED RESEARCH 


• KAMAN AEROSPACE CORPORATION - AUTOMATED MATH MODEL 
IMPROVED FOR MATCHING EXPERIMENTAL DATA. 


I INSTITUTE FOR COMPUTER APPLICATIONS IN SCIENCE AND 
ENGINEERING - IDENTIFICATION OF EQUIVALENT PDE SYSTEMS 
TO MATCH MEASURED DATA. 


SUMMARY 


• COMPUTER PROGRAM FOR REDESIGNING STRUCTURAL MODES TO 
REDUCE RESPONSE HAS BEEN INITIATED. 

• LINEAR REGULATOR APPROACH IN MODAL COORDINATES HAS BEEN 
IMPLEMENTED. TRANSFORMATION OF SOLUTION TO PHYSICAL 
STRUCTURE IS A MAJOR PROBLEM. 

• SOLUTION OF STIFFNESS EQUATIONS AND DAMPING EQUATIONS 
CAN BE DONE SEPARATELY AS NXN SET OF (MATRIX RICCATI) 
EQUATIONS. 

PLANNED EFFORT FOR '82 

• INCLUDE MASS OF CONTROL 

• STUDY WEIGHTING TO MINIMIZE OR SELECT CROSS-TERMS 
I IMPLEMENT PHYSICAL COORDINATE SOLUTION 

• STUDY POTENTIAL FOR "BENEFICIAL" CROSS TERMS 


149 





ALGORITHMS FOR ON-LINE 


PARAMETER AND MODE SHAPE ESTIMATION 


Frederick E. Thau 
The City University of New York 
New York, N.Y. 


One of the schemes that has been proposed for the adaptive control of la.rge 
flexible space structures la shown In figure 1. This approach Is based upon a 
modal deconiposltlon of the dynamic response of the flexible structure and Is 
designed to make uae of the parallel processing features of modern minicomputers. 
Satisfactory performance of the parallel structure Identification technique shown 
In figure 1 can he achieved only when the approximation functions noted In the 
figure correspond to the natural modes of the flexible structure. The work 
su mma rized here presents a technique for estimating both mode shapes and modal 
parameters. 


m 

w ( s, t ) = q. ( t, p. ) c j ( s ) + V”( s, t ) 

i=l ' ‘ ' 

(a) Motion of flexible structure. 


q(k + l) = A,q (k) + A^qU-l) + Bj U ( k ) + B2U(k-l) 


y(k) = (k» + w(k) 

Q 

(b) Modal description. 


Figure 1.- Problem formulation. 
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Figure 2 shows the analytic background for the formulation of the on-line 
' identification problem. The jnotion of the flexible structure is expressed as a 
sum of NM teinns involving the :modal shape function E^Cs) and the modal amplitudes 
qiCtjPi^I where s denotes the spatial variable and t represents time. The 
parameter vector pi represents a set of parameters l^^±y^2±9^1±y^2±^ each 

mode. The identification problem is formulated in terms of the modal description 
shown in figure 2(B), where the Ai,i=l,2 are diagonal matrices of order NM, U(k) 
denotes a vector of actuator signals and the are rectangular matrices. Matrix 

H, relating the modal anplitudes to the measurement vector y(k), has columns that 
are linearly independent in the modal model. On-line measurements {y(k)},{U(k)} 
are to be processed in order to obtain estimates of the matrices A^, and H. 




Figure 2.- Distributed adaptive control. 
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Algorithms for on-line parameter identification are shown in figure .3. 
Assuming that the approximation functions are given, the modal parameters are 
updated using an output error formulation where the error e^ is given by 


e^(k-l) = q.(k-l) - A^^q^(k-2) - A2^q^(k-3) - Si^F^(k-2) - B2iF^(k-3) 


where Fi(k) 
satisfy 


denotes a modal force coiiq>onent. 


The weights 


W. 


are selected to 


Wiq^(k-2) + W2qJ(k-3) + W3Fj(k-2) + W^F^(k-3) < 2 


to insure stability of the identification algorithm. 

The approach for updating the approximation functions is based on combining 
the modal equations shown in figure 2 into the single relation involving matrices 
and where 






i = 1,2 


and (H)^ denotes the pseudo-inverse of H. Using a regression analysis approach, 
data are collected over a time interval of length N and a least-squares estimate 
of Mi and Ni is obtained. When the number of modes in the modal approximation 
is the same as the number of sensors, the eigenvectors of Mi correspond to the 
columns of H. 


Pj(U 


'W,q. (k-2)' 

p. (k-f) + e. <k-f)- W^q. (k-3) 

W^F. (k-2) 

W.F. (k-3) 
4 I j 


(a) Pole-zero characteristics - output error formulation. 


y(k + I) = M,y(k)+M 2 y(k-I) + N|U(k) + N 2 U(k-ll + n(k + l) 


Ntv(N) = S(NIM^ + V(N): = 


( N ) Y ( N ) 


(b) Approximation functions - regression analysis. 


Figure 3.- On-line parameter identification. 
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A schematic diagram of the flexible beam used to test the Identification 
algorithms is shown in figure 4. Sensor and actuator locations along the beam are 
shown along with the modal frequencies and mode shapes for the first three flexible 
modes obtained from the SPAR analysis program. 


SENSOR 


ACTUATOR 


1.81 H 


4.99 Hz 


9.74 Hz 


Longitudinal S tat ion , f t , 


Figure 4.- Modes of interest. 







Figure 5 shows the modal decomposition for mode 3 unforced response. A set 
of approximation functions, obtained from the SPAR analysis was used and eight 
modes, were assumed in the real-time program that produced the modal amplitude 
signals shown. For this response there is some excitation of modes 1 and 2. 
However, higher frequency modes do not appear to be excited. 




r i I i J 

0 10 20 30 40 


Tlme» sec. 


Figure 5.- Modal decomposition. 
Eight SPAR modes; mode 3 excited. 



Outputs> of four on-line ident 
of the third vibration mode are sh 
oacillation occurs, the identifiers 
initiation of the susiained oscill 
rapidly identified. The parameter 
the matrix Ai, i = 1,2, 
modal parameters derived 



Figure 6.“ Output of par 







figure 7 show.a the outputs of four on-line identifiers resulting from a 5-Hz 
excitation produced by actuator 1; Note that the parameters of the first vibration 
mode are Identified but the high, frequency modal parameters are not found as a 
result of the 5-Hz excitation. 


mil 




nil 

wiirr 


Time, sec 


Time, sec 


Figure 7.- Output of parameter identifiers. SPAR derived values 
Sinusoidal excitation by actuator 1 at 5 Hz. 







j Figure 8 demonstrates the result of applying the mode shape identification 

procedure (a) to a simulation wherein the beam was given an arbitrary displacement 
and (b) in an experiment wherein only the first flexible mode was excited. In the 
simulation the normalized mode shape obtained agreed with that produced by the 
SPAR analysis. However, in the laboratory experiment a distortion in the identified 
mode shape was observed. In both cases the mode shape was identified on line after 
approximately 1 second of data processing. 



(b) Experimental 



Longitudinal station, ft 


Figure 8.- Identification of mode shape. 
First flexible mode. 
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A table of identified modal parameter values that resulted from a simulation 
incorporating both the mode shape and modal parameter identification schemes is 
shown in figure 9. Identification of mode shapes and parameters was accomplished 
in approximately 1 second of processing measurements from eight sensors distributed 
along the beam. 


[simulation - all modes excited at parameter values] 


Flexible mode 

Initial 

Final 

Actual 

1 

-0.8495 

1.8726 

1.87267 

2 

-1.7194 

1.0991 

1.099181 

3 

0.09918 

-0.71936 

-0.71940 

4 

0.83907 

-1.9960 

-1.9960 

5 

0.8726 

0.15047 

0.15049 

6 

-2.996 

1.8387 

1.83907 


Figure 9.- Identification of modal parameters. 


SUMMARY 


• Presented algorithms for on-line parameter and mode-shape estimation 

• Examined identification performance using computer simulations and a 

limited number of laboratory experiments 

• Future work will include: 

(1) Further experimental studies 

(2) Development of separable nonlinear least-squares approach 

to identification 

(3) Development of on-line performance measures for sequential 

processing decisions 
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COMPONENT NUMBER AND PLACEMENT 


IN LARGE SPACE STRUCTURE CONTROL 


Wallace E. VanderVelde 
Department of Aeronautics and Astronautics 
Massachusetts Institute of Technology 
Cambridge, ^lassachusetts 



PURPOSE: 


Provide an objective means of assisting the designer 

OF THE CONTROL SYSTEM FOR A LARGE FLEXIBLE SPACE 
STRUCTURE IN HIS CHOICE OF HOW MANY ACTUATORS AND 
SENSORS TO INCORPORATE IN THE SYSTEM^ AND WHERE TO 
LOCATE THEM ON THE STRUCTURE. 


WHATHENEEDIS ; 

0 A QUANTITATIVE MEASURE OF HOW WELL A SYSTEM CAN BE CONTROLLED 
WITH A SPECIFIED SET OF ACTUATORS 

0 A QUANTITATIVE MEASURE OF HOW WELL A SYSTEM CAN BE OBSERVED 
WITH A SPECIFIED SET OF SENSORS 

0 A MEANS OF RECOGNIZING THE EFFECTS OF COMPONENT FAILURES IN 
THESE MEASURES 

0 A MEANS OF OPTIMIZING THE LOCATIONS OF ACTUATORS AND SENSORS 
SO AS TO MAXIMIZE THESE MEASURES 
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A MEASURE OF THE DEAREE OF CONTROLLABILITY 


1. Find the minimum control energy strategy for driving the system 

FROM A GIVEN INITIAL STATE TO THE ORIGIN IN A PRESCRIBED TIME. 

2. Define the region of initial states which can be returned to the 

ORIGIN WITHIN SPECIFIED LIMITS ON CONTROL ENERGY AND TIME USING 
THE OPTIMAL STRATEGY. 

3. Define the Degree of Controllability to be some measure of the 
SIZE OF THIS region, 


STEP 1. MINIMUM ENEROY CONTROL 


1 I T 

Problem statement: Min ^ ~ ^ I H. Rudt 

o 

subject to X “ ^21 

x(0) , T given 

x(T) = 0 

Solution: u(t) = R*"^ 0pp ^xp 

The (t) are partitions of the transition matrix for the system 


_d^ 

dt 



-1 T 

-BR 


-A 
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STEP 2 . 


THE RECOVE R-Y REfilON 


Using optimal control, the cost to return to the origin can be 
expressed : 

E = i x(0)^ V(0) x(0) 

with V = AV + VA*^ - BR~^ b'^ V(T) =0 


For specified values of time, T, and control energy, , the recovery 
region is the interior of the space bounded by the surface 


x(0)^ V(0)“^ x(0) = 2 Eg 


STEP 3. THE SIZE OF THE RECOVERY REF ION 


First scale the state variables such that equal displacements in all 
directions are equally important, 


2 = Dx 





where x. is the minimum initial value of x. one would like to be able 
^min ^ 

to drive to the origin with constrained time and control energy. 
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STEP 3 (continued) 


Then define a weighted measure of the volume of the recovery region 
in the scaled space. 




Vg = volume of the largest ?jphere which can be inscribed in the 
elliptical boundary of the recovery region 

Vp = volume of the recovery region 



V - V 
R S 


weighted volume. 


Also 


DC= 

v„ ^ (A, ) 


TT ^ -1/2 


where the are the eigenvalues of V(0)~^ D 


Alternatively 


^ TTi (l/i) 


1/2 


where the are the eigenvalues of DV(0)D. 


Then 




n/2 


An analytic solution is available for V(0) , 
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DEGREE OF CONTROLLABILITY 
FOR n FREE -FREE BEAM 


S First flexible mode only 



DEGREE OF CONTROLLABILITY 
FOR A FREE-FREE BEAM 


Second flexible mode only 




DECREE OF CONTROLLABILITY 
FOR fl FREE-FREE BERM 


First and second flexible modes 



DEGREE OF CONTROLLABILITY 
FOR A FREE-FREE BERM 


First and second flexible modes 




DEGREE ’OF CONTROLLABILITY 
FOR A FREE-FREE BEAM 


rsi 



ONE MORE CONSIDERATION 


We want the measure of controllability to reflect the fact 

THAT A SYSTEM WITH MORE ACTUATORS OF EQUAL EFFECTIVENESS HAS 
GREATER CONTROL CAPABILITY THAN ONE WITH FEWER ACTUATORS. 

The DC just defined is made proportional to the number of 

ACTUATORS PLACED AT THE SAME LOCATIONS IF THE ELEMENTS OF R 
ARE SCALED INVERSELY WITH THE NUMBER OF ACTUATORS. 

For DIAGONAL R^ CHOOSE Roii TO REFLECT THE RELATIVE COST OF 
THE DIFFERENT CONTROLS. ThEM 


m = total number of actuators 



A_ A S U RE OF THE DECREE OF 
OBSERVABILITY 


1. In terms of an information matrix^ determine how much 

INFORMATION CAN BE DERIVED ABOUT THE SYSTEM STATE IN 
TIME STARTING FROM ZERO INFORMATION^ USING THE GIVEN 
SET OF SENSORS. 

2, Define the Degree of Observability to be a measure of 

THE size of this INFORMATION MATRIX. 


STEP 1. THE INFORMATION MATRIX 


We want the Degree of Observability to be a property of the system^ 

NOT OF THE ENVIRONMENT IN WHICH IT OPERATES. So DO NOT CONSIDER 
STATE DRIVING NOISE, 

Then 

* T T -1 

J = “ JA - A J + C M -^C 

J(0) = 0 

CORRESPONDING TO THE SYSTEM MODEL 

X = Ax + Bu 
^ = Cx + n 

= N(5(t2“tj^) 

An analytic solution is available for J(T). 
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STEP 2. THE SIZE OF THE MATRIX 


One way to measure the size of J(T) is to indicate the volume 

CONTAINED WITHIN THE SURFACE 
J(T)"^v = 1 

But the variables should be scaled to reflect the relative importance 
OF errors in the different state variables, 


w = F V 



WHERE e. IS THE MAXIMUM TOLERABLE ERROR IN THE ESTIMATE OF 
max 


The Degree of Observability is defined with respect to this volume 
IN the space of equally important errors (w) just as the Degree of 
Controllability was defined for the volume of the recovery region 

IN THE SPACE OF EQUALLY IMPORTANT CONTROL CHARACTERISTICS, 

^ '''r - 

''s = Vi 

WHERE THE l^i ARE THE EIGENVALUES OF FJ(T)F. 
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DEGREE OF OBSERVABILITY 
„ FOR A FREE-FREE BERM 

o 



DEGREE OF OBSERVABILITY 
FOR A FREE-FREE BERM 

First and second flexible modes 




RECOGNITION OF C0^1PnNENT FAILURFS IN 
THESE MEASURES 

Let f INDICATE THE STATE OF FAILURES AMONG THE ACTUATORS OR SENSORS. ThEN 
FOR EVERY DC(f) OR DO(f) CAN BE COMPUTED AS JUST DESCRIBED. 

From the statistical model for failures of the different components one can express 

P[f(t) = f^]. 

The average^ over the mission period^ of the expected Degree of Controllability or 
Observability is taken as the final measure. 

T 

1 r ^ _ 

AD C = ^ / DC(t) dt 
"■o 

= i f " E DC(f.) 

m J i 
o 

- E DC(f,) if""' 

i mJ 

o 

The Degree of Observability is computed in the same way. 


First and second flexible modes 
2 actuators: 1 FIXED AT END 



P[f(t) = f.] dt 
P[f(t) = dt 
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OPTIMIZING C 0 M P 0 N E N T 1 0 C A T I Q N S 

The average Degree of Controllability is a function of the choice of 

ACTUATOR LOCATIONSJ LET INDICATE THE VARIOUS ADMISSIB LE CHOICES, 

ADC ^ ADC (^) 

An optimization routine is required to find the OPTIMUM LOCATIONS. 

ADC* = max ADC (jE,) 

The optimum sensor locations are defined in the same way using the 
AVERAGE Degree of Observability. 


First and second flexible modes 
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CHOICE OF THE NUMBER OF COMPONENTS 


The optimum average Degree of Controllability is a function of 
THE number of ACTUATORS IN THE SYSTEM, 

ADC* ^ ADC* (m) 

With the likely constraint that multiple actuators cannot be placed 
IN the same location^ ADC*(m) SHOWS DIMINISHING RETURNS WITH 
increasing m. 

The same is true for the optimum average Degree of Observability. 
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A PROGRAM PLAN FOR THE DEVELOPMENT OF 


FAULT TOLERANT LARGE SPACE SYSTEI4S 


Paul Motyka 

Charles Stark Draper Laboratory, Inc* 
Cambridge, Massachusetts 


Objectives 


for fault tolerance In LSS 

lue characteristics of LSS which 
lerance 

:atus of fault tolerant systems 


im plan to validate and demonstrate 
cult tolerance for LSS 


e Need for Fault Tolerant LSS 



Number of System 
Comnonents 


N « 200 





Characteristics of LSS That Affect Fault Tolerance 


• Dimensionality 

- Will have to consider a much larger number of components 
than In any previous application 

- May affect the achievement of real time operation and 
computational accuracy 

• Precise accuracy and stability requirements 

- Small failures will have to be detected and Isolated quickly 

• Structural mode and physical displacement effects on sensors 

- May be comparable in magnitude to the failures which 
must be detected 

• Environmental effects 

- Dynamic effects such as large angle slewing maneuvers 

- Thermal effects may produce chanoes In modal 
characteristics 

• Spillover (Model Order Reduction) effects 

- May Introduce uncertainties which can be falsely 
Interpreted as failures 

• Many diverse types of sensors present 

- Must be collectively accounted for if system wide 
fault tolerant capability Is to be achieved 

• Multiple experiments operating on a single LSS 

- Interactions among experiments may result In false alarms 

- A robust fault tolerant system may be required to tolerate 
changes In modal characteristics with experiments 

• The modal frequencies may be In the controller bandwidth 

- Filtering to reduce modal effects and Improve fault 
tolerant capability mav not be possible 
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Status of Fault-Tolerant Technology for LSS 


Little work has been done In this area 

A large body of knowledge concerning fault tolerant systems 
has evolved In spacecraft and avionics applications 

This material forms o solid foundation for the development of 
fault tolerant technology for LSS 


Elements of a Program Plan to Validate 
the Concept of Fault Tolerance for LSS 

System modeling 

- Generate an analytic model of the LSS 
Environment modeling 

- Define LSS tasks^ maneuvers and disturbances 
Requirements definition 

- Define LSS accuracy, stability and relloblllty 
requirements 

Component modeling 

- Define anolytlc models, error effects, noise, flexibility 
effects, etc. 

- Uncertainties establish the fault tolerant ccpablllty of LSS 
Fault tolerant techniques development 

- Develop algorithms to detect and Isolate faults and 
reconfigure LSS 

Simulation development 

- Both nonreal time and real time capabilities will be 
needed 

Development of o foult tolerant data processing copoblllty 

- Needed to carry out computations associated with the LSS 

- A firm basis for this technology exists 



Elements of a Program Plan to Validate 
the Concept of Fault Tolerance for LSS 
(concluded) 


Conceptual design 

- Preliminary definition of algorithms^ components/ 
architectures/ etc. and assessment of design 
alternatives 

System design 

- Detailed and specific determination of parometers/ 
components/ architectures and algorithms for selected LSS 

Preliminary Implementation 

- Definition of system software and hardware 

- Partitioning of functions amona subsystems 

Validation of fault tolerance concept and demonstration 
of capability 

- Use simulation of complete LSS 
Requirements evaluation 

- Use simulation to assess ability of fault tolerant 
LSS to meet accuracy and stability requirements 

- Analytic techniques must be employed to evaluate 
the reliability of the LSS 


Conclusions 


There Is a definite need for fault tolerance In LSS 

LSS have unique characteristics which Impact fault tolerance 

Very little work has been done regarding fault tolerance for 
LSS although a solid base exists from spacecraft and avionics 
applications 

A program plan for the validation and demonstration of the 
concept of fault tolerance for LSS has been developed 




LARGE SPACE STRUCTURES 
CONTROL ALGORITHM CHARACTERIZATION 


E. Fogel 

Charles Stark Draper Laboratory, Inc. 
Cambridge, Massachusetts 


Computation Consideration in 
LSS Control/Identificotion 

• Algorithms 

• Structures 

• Computation considerations 


Spillover Effect 
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MODEL 


X = AX + BU 


y = CX 


A = diag Hj 



1 

0 


Separation to : 


controlled modes 


: residual modes 


X = 




y = C(.X^. + CfjXp 


B 


c 


X + 



u 


LAC/HAC 

LAC ; local feedback colocated sensor/octuotor pairs 

Augment damping 

HAC ; dynamic feedback to control a reduced order model 

* frequency shaped K.F. 
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HAC/LAC Control Algorithm 

LA£; Ul = G y G = ? > 0 

MAC ; = Q f + M 

\ = -+ BU + K(y - CjJXj.) 

Uh = GiXc + 62/' 

U = Ul + 

Rate HAC rote = 1/2 LAC rote 


LAC/HAC BASED COMPUTATION REQUIREMENTS 



I I 

COMPUTER 
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LOG APPROACH 


Solution 


solve LQ6 for 

J - /[lUc IIq^ + ||u ||„2] dt 


Implementation 
U = T K 

% - AcXc + BcU+KT(y-CcXc) 

K » - R’* bJp 

PAj. + aJp + Q - ptBcR'^tbJp = 0 


T; BrT = 0 

f : T Cr - 0 


Bj.T i 0 orthogonality 

TC- * 0 conditions 

I# 


Closed loop: 



Xc 


. — 
Ac - BJK 0 BJK 



d 

^R 

e 

- 

0 Af^ 0 

0 0 Ac - SCc 


Xr 

e 



185 












Beam Instrumentation 

5 clusters of x y f accelerometer (sensor) 

L proof mass (actuator) 
at: TOP/ Middle/ Bottom 

DIMENSION (m) of SENSOR VECTOR' AND ACTUATOR VECTOR = 3 X 2 » 6 


187 



60m REFLECTOR 



V *■ 


OFF SET FEED 
F/O - 2.0 

NONMETALLiC FEEDMAST 
TAPERED TENSION TRUSS 



I 

50m 


Antenno Instrumentation (detail) 


13 CLUSTERS OF COLOCATED SENSORS/ACTUATORS AS FOLLOWS: 
1. HAST/ORBITER ATTACHMENT: 


SENSORS 2 DOF ACCELEROMETER PKG (x, y) 

ACTUATOR 2 DOF PROOF MASS 

2. REFLECTOR HUB (WHERE FEED SUPPORT MAST IS ATTACHED TO ANTENNA SUPPORT MAST) 


SENSORS 2 DOF 

1 DOF 

ACTUATOR 2 DOF 

1 DOF 


ACCELEROMETER PKG (x, y) 
RATE GYRO (TORSION AXIS) 
PROOF MASS PKG 
TORQUE WHEEL 


5. 8 CLUSTERS OF INSTRUMENTS AROUND RIM OF REFLECTOR: 

SENSORS 2 DOF ACCELEROMETER (TANGENTIAL,+ z) TENSIOMETERS ON GUY 

WIRES 

ACTUATORS 2 DOF PROOF MASS ( tangent lol +Z), Guy Tensioner 

g. MIDDLE OF FEED SUPPORT; 

SENSORS 2 DOF ACCELEROMETER PKG (x^y) 

ACTUATORS 2 DOF PROOF MASS (X,y) 



Antenna Instrumentation (cent.) 


5. FEED HAST/SUPPORT MAST ATTACHMENT: 


SENSORS: 

2 DOF 

ACCELEROMETER 

(x,y) 


1 DOF 

RATE GYRO 

(TORSION) 

ACTUATORS 

2 DOF 

PROOF MASS 

(x,y) 


1 DOF 

TORQUE WHEEL 

(TORSION) 

AT FEED: 

SENSORS 

2 DOF 

ACCELEROMETER 

(y,z) 

ACTUATORS 

2 DOF 

PROOF MASS 

(y,z) 


• DIMENSION OF SENSOR/ACTUATOR VECTORS (M) = 2+3+2+3+2+8x3 = 36 


LQG A ND HAC/LAC COrPUTATlOWAL SIZING 

• THESE ALGORITHMS HAVE BEEN SIZED IN TERMS OF 

• FLOATING POINT OPERATION (FLOP) DEMANDS 

• STORAGE FOR VARIABLES 

• INPUT/OUTPUT DATA FLOW 

• FLOP SIZING (PER CONTROL CYCLE) DONE AS A FUNCTION 
OF THE NUMBER OF CONTROL STATES AND THE NUMBER OF 
SENSOR/ACTUATOR PAIRS 

• STORAGE FOR VARIABLES AND I/O SIZING DONE FOR 
SPECIFIC STRUCTURE EXAMPLES 


Input/Output Data Flow Rates 


Assumption • Control bandwidth 50 Hz 
• Accuracy - 2 byte/word 


[ 


Sampling frequency 
Command frequency 


250 Hz 


Data Flow rate 


per 


sensor 

actuator 


500 [Bytes/seg 


Total : Beam: 

Antenna: 

1553B bus capacity 


3^000 tJytes/sec] 
18,000 Uytes/sec] 

A8,000 [Bytes/secl 
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LQ£JiZINfi 



BEAM 

MIENM 

SENSOR/ACTUATOR PAIRS On) 

6 

36 

CONTROL STATES (n^) 

20 

20 

FLOP PER CYCLE* 

1A20 

AA20 

VARIABLES** 

752 

2312 

I/O PER CYCLE 

12 

72 


*INCLUDES SENSOR COMPENSATION FLOP (120 FOR BEAM. 720 FOR ANTENNA) 
**INCLUDES SENSOR COMPENSATION VARIABLES (60 FOR BEAM. 360 FOR ANTENNA) 


FLOP/CYCLE 

LOG ALGORITHM: 2n^+n +4n m 




\i\ 





SENSOR/ACTUATOR PAIRS (m) 

6 

36 

CONTROL STATES (n^) 

12 

12 

FLOP PER CYCLE* 

633 

A608 

VARIABLES** 

570 

3060 

I/O 

12 

72 


INCLUDES SENSOR COMPENSATION FLOP (120 FOR BEAM, 720 FOR ANTENNA) 
INCLUDES SENSOR COMPENSATION VARIABLES (60 FOR BEAM. 360 FOR ANTENNA) 


FLOP/LAC CYCLE 

(LAC FLOP + H HAC FLOP): + 5**™* 


TYPICAL AP 
CAPACITY e 
250 LAC/s 


STATES UNDER HIGH-AUTHORITY CONTROL (n^) 








COHPUTATIOW LOAD 


Structure 

Algorithm 

m 

n 


£a£fi 

K Flops/sec 

% 6PC 
capacity 

X typical 
AP capacity 

Beam 

LQG 

6 

12 

6 

55 

67 

7 



6 

16 

10 

112 

135 

14 


HAC/LAC 

6 

12 

6 

3A. 

<ll 

4 



6 

16 

10 

57 

69 

7 

Antenna 

LQG 

36 

A2 

6 

260 

300 

36 


HAC/LAC 

36 

^2 

6 

750 

900 



in = # of sensors/octuotors 
n - # of inodes In model 



# of controlled inodes 


SYSTEH-IJEN I IFICATXQtl. COMPUTATIONAL SIZING 

• ARMA-LEAST SQUARES ALGORITHM SIZED FOR FLOP AS FUNCTION OF 
MODEL ORDER (n) AND NUMBER OF SENSOR/ACTUATOR PAIRS (m) 

• FLOP REQUIREMENTS FOR THIS ALGORITHM ARE SO LARGE THAT 
IMPLEMENTATION IN A FLIGHT SYSTEM OR ITS GTF ANALOG IS 
PRECLUDED 

• EVEN IMPLEMENTATION IN GROUND-BASED COMPUTERS IS 
CONSIDERED QUESTIONABLE. BUT THIS STUDY ASSUMES 
A GROUND-BASED IMPLEMENTATION 


• NOTE: SOME OTHER SYSTEM IDENTIFICATION ALGORITHM MAY BE 
IMPLEMENTABLE IN A FLIGHT SYSTEM 
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Alaorlthm Assessed - Least Squares 
Motivation for chooslna LS 




• Relative high spectral resolution 

• Comparable to other alnorlthm In computation 
complexity 

e.g,: Covariance algorithm 

: Maximum Entropy 

• "Better" algorithms - considerably more complicated 

• Less complex algorithms - considerable penalty In 

performance 

• LS - robust to order reduction 

• Useful for - control design 

- self tuning regulators 


Identification Algorithm Sizing 

Assume the ARMA model 

^'k “ ^ ^l^k-l ^ ®lVl 

where y|^ = vector of measurements (sensors) 
at cycle k 

U|^ = vector of control Influence at cycle k 
we can write 

^k [^1 • ;^n 1 ®o i ■ ’ • ®n-lj 


= a.Zk 

Use least squares Identification 


I 

yk-n 

Uk 

I 

^k-n+1 


S YS TErUMTIFICAUOfl. ALGORITHM. FLOP R EQUIREWENTS 


SENSOR/ACTUATOR PAIRS (m) 
MODES MODELLED (n) 

OFF-LINE MEGAFLOP 
FOR «lOOO CYCLES 

OFF-LINE FLOP/CYCLE 

OFF-LINE MEGAFLOPS 
(a 250 CPS) 

ON-LINE FLOP/CYCLE 

ON-LINE MEGAFLOPS 
(8 250 CPS) 


BEAM 

ANTENNA 

6 

36 

12 

92 

35A.2 

297,779 

88,552 

79,999,881 

22.1 

18,611 

169, 78A 

73,601,179 

92. 5 

18,900 


.01 


AVIONICS DATA PROCESSING 
throunhput (MFLOPS) 
0.1 1-0 


10.0 


100 


Existing Space Processor 


Near-Term Processor 


C/l o z o . 

O e» t/> Tl -D I 

-o O </) I I XT 

I I O *— U1 

INJ xr I N) O 

M o> t 


P -Processors 


Projected Trends 
(avionics) 


Existing Control 

LSS Control (ATB Models) 

LSS Identlflcotlon 
Other Space ROTS 






ro 


UD 

i£> 

o 

o 


3 5 s 

0 3> -n 

O I 

1 N) 

4=- ro 

4=“ 

3) 


V/1 


oo 

o 


oo 

<J 1 


i£> 

o 




(D 


[Wocjel Dependent j 


Housekeeping 


Payload Processing y ,Slgnol Process Ing^ 
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PARTITIONING OF LARGE SPACE STRUCTURES 


VIBRATION CONTROL COMPUTATIONS 


J . Kernan 

Charles Stark Draper Laboratory, Inc 
Cambridge, Massachusetts 


LQG SIZING 



BEAM 

ANTENNA 

SENSOR/ACTUATOR PAIRS (m) 

6 

36 

CONTROL STATES (n^,) 

20 

20 

FLOP PER CYCLE* 

1A20 

mo 

VARIABLES** 

752 

2312 

I/O PER CYCLE 

12 

72 


'includes sensor compensation flop (120 FOR BEAM, 720 FOR ANTENNA) 
"includes sensor COMPENSATION VARIABLES (60 FOR BEAM, 360 FOR ANTENNA) 


HAC/LAC SIZING 



BEAM 

ANIENNA 

SENSOR/ACTUATOR PAIRS (m) 

6 

36 

CONTROL STATES (n^.) 

12 

12 

FLOP PER CYCLE* 

633 

A608 

variables" 

570 

3060 

I/O 

12 

72 

INCLUDES SENSOR COMPENSATION FLOP 

(120 FOR BEAM. 

720 FOR ANTENNA) 


INCLUDES SENSOR COMPENSATION VARIABLES (60 FOR BEAM, 360 FOR ANTENNA) 
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niSTRlRIITinN OF V lRRATinN r.ONTRQL AND SENSOR COMPENSATION COMPUTATIONS 


• A RANGE OF CHOICES BETWEEN THE TWO FOLLOWING 
EXTREMES WAS INVESTIGATED; 

CENTRAL COMPUTATION OF BOTH CONTROL 
AND SENSOR COMPENSATION 
• CONTROL COMPUTATIONS DISTRIBUTED AMONG 
CENTRAL AND LOCAL PROCESSORS; LOCAL 
PROCESSORS ALSO PERFORM SENSOR COMPENSATION 
• PORTIONS OF THE CONTROL COMPUTATIONS CAN 
BE DISTRIBUTED BECAUSE THE ROWS OR COLUMNS 
OF THE MATRICES INVOLVED CORRESPOND TO 
INDIVIDUAL SENSORS OR ACTUATORS 
• DISTRIBUTION IS SUPPORTED BY THE FACT THAT 
THE MATRICES INVOLVED ARE EITHER CONSTANT 
OR INFREQUENTLY CHANGED 


LSSC COMPUTATION DISTRIBUTED AMONG CENTRAL AND LOCAL PROCESSORS 


9999 ^ 


-SENSOR 

CLUSTER 


RAW SENSOR READINGS ^ 


COMPENSATION FOR SCALE 
FACTOR, BIAS, AND DRIFT 

COMPENSATED READiNG^ 


— 


"v" 



- 

X X X X 


Y 



Y' 

X X X X 


V 


MXV">* 

Y' 

X X X X 


V 



Y' 

— _ 



— 


KALMAN GAIN COLUMNS 
FOR THESE 4 SENSORS 


WEIGHTED 

SENSOR 

READINGS 





rvu 


“y'" 

LyJ^ 

MX^>— ► 

Y' 

V' 



¥ 


CENTRAL 

PROCESSOR 

ACTIVITIES 


VIBRATION 

CONTROL 

ALGORITHM 


4-ACTUATOR 

CLUSTER 


9999 


ACTUATOR COMMANDS 


FEEDBACK GAIN ROWS 
FOR THESE 4 ACTUATORS , 


XXX 

XXX 

XXX 

XXX 


\ 







MXV 


LOCAL PROCESSOR OUTPUT ACTIVITIES 



ESTIMATED | 
STATE 
VECTOR 


ETC. 


ETC. 
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I^COimDATl.Q«S-F.QR..CENTl^ 


LOG COMPUTATIONS THAT SHOULD BE DISTRIBUTED TO LOCAL PROCESSORS^ 

IN ORDER OF DECREASING PREFERENCE 

. COMPENSATING SENSOR READINGS FOR SCALE FACTOR^ BIAS/ AND DRIFT 
. APPLYING FEEDBACK GAIN TO ACTUATOR COMMANDS 
. APPLYING KALMAN GAIN TO COMPENSATED SENSOR READINGS 

HAC/LAC computations that should be DISTRIBUTED TO LOCAL PROCESSORS/ 

IN ORDER OF DECREASING PREFERENCE 

. COMPENSATING SENSOR READINGS FOR SCALE FACTOR/ BIAS/ AND DRIFT 
, APPLYING LOW-AUTHORITY GAIN TO ACTUATOR COMMANDS 

, APPLYING HAC FEEDBACK & FILTER GAINS TO ACTUATOR COMMANDS 

. APPLYING KALMAN GAIN TO COMPENSATED SENSOR READINGS 

FREQUENCY SHAPING FILTER COMPUTATION (HAC/LAC) SHOULD STAY IN 
CENTRAL PROCESSOR 


LOG ALGORITHM (n^ - 20) FOR BEAM 




I 


CENTRALIZED 


DISTRIBUTED 




LOG ALGORITHM (n^ - 20) FOR ANTENNA 




CENTRALIZED 


DISTRIBUTED 


HAC/LAC ALGORITHM (n^ « 12) FOR BEAM 




•THESE OUTPUTS ARE SENT ON 
EVERY OTHER CYCLE 


CENTRALIZED 


DISTRIBUTED 
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HAC/LAC ALGORITHM (n^ - 12) FOR ANTENNA 



•THESE OUTPUTS ARE SENT 
ON EVERY OTHER CYCLE 



CENTRALIZED DISTRIBUTED 

CONCLUSION 

• VIBRATION CONTROL OF LARGE SPACE STRUCTURES 
IS COMPUTATIONALLY DEMANDING - DRIVEN BY 

. NUMBER OF VIBRATION MODES CONTROLLED 

. NUMBER OF SENSOR/ACTUATORS PAIRS 

. CONTROL BANDWIDTH 

• DISTRIBUTION OF THE VIBRATION CONTROL COMPUTATIONS 
AMONG CENTRAL AND LOCAL PROCESSORS CAN SIGNIFICANTLY 
REDUCE THE THROUGHPUT REQUIRED FROM THE CENTRAL 
PROCESSOR AND MAY ALSO RESULT IN IMPROVED PERFORMANCE 
DUE TO REDUCED TRANSPORT LAG 
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SHUTTLE FLIGHT CONTROL 
AND STRUCTURE INTERACTION 


Michael Paluszek 

Charles Stark Draper Laboratory, Inc 
Cambridge, Massachusetts 


• A brief overview of the DAP 

• Results of the dynamic Interaction study of 
the Orblter with the General Dynamics Beam 
attached 

• Preliminary results of the three degree of 
freedom payload parametric study 


This is a diagram of the Orbiter Flight Control System (FCS) . The functions 
relating only to the QMS have been eliminated for clarity. The vehicle can be 
maneuvered either manually or automatically. To perform an automatic maneuver the 
pilot specifies a new orientation for the Orbiter. This is processed by the steering 
processor, which sends a corresponding attitude error to the RCS processor. Jets 
are commanded to fire to reduce this error. The jet firings produce attitude changes 
that are measured by the IMU and sent to the state estimator at 6.25 Hz. The state 
estimator produces rate and acceleration estimates, filtering out high-frequency 
oscillations. The steering processor produces a new rate and attitude error from 
this data. 



I : 1 


FCS 
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Characteristics of the Orblter and DAP 
Orblter 

• Low Inertial cross coupling 

• First bending frequency 0.431 Hz 

• Large RCS Jet coupling particularly roll to yaw 

• - .883x10® slug-ft^, 12. 5X of 

• Operational maneuver rate limit of 2 °/sec 
DAP 

• Phase Plane autopilot assuming decoupled axis 

• State Estimator uses IMU data only 

• second order filter characteristic 

• 6 dB down at .06 Hz with vernier gains 

• Designed assuming a rigid vehicle 

This flow chart maps the interrelation of different criteria for the Orbiter 
flight envelope. The three parameters that are used to define an envelope are 
payload pointing accuracy, propellant budget and payload loads. While often con- 
sidered independently, they affect one another as indicated in the block diagram. 


Definition of the Orblter Flight Envelope 
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Classes of System Response 
for the Orblter with o Flexible Payload 

1. The Autopilot does not respond to the payload oscillations 

• vehicle oscillations small 

• paylood oscillations may be large 

2. The Autopilot responds to payload oscillations 

• vehicle rate or attitude errors exceed set limits 

• vehicle oscillations may diverge 

• attitude excursions may be small 


This graph defines the regions of autopilot interaction for the Orbiter with a 
flexible payload. The x axis is the roll fundamental frequency and the y axis is the 
ratio of payload roll Inertia to empty Orbiter roll inertia. Little or no interaction 
occurs in the region to the right of the solid line. In the region to the left some 
closed-loop response has been observed. The dotted line indicates the state estimator 
6 dB point. The ruled region is the current area of interest for the parametric 
study of the DAP response. 
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This table gives the results of the dynamic interaction analysis of the 
Orbiter/GS beam combination* The columns list the maneuver performed, the fuel 
expended, the type of Orbiter and payload oscillations, whether the DAP sent 
oscillatory firing commands to the jets, and comments. 

Except for the 5® roll case with rate limit clampdown, the vehicle and pay- 
load remained stable. In cases 4 and 7 through 10 oscillatory commands were sent 
to the jets indicating vibration feedback. Note that 4 was the worst case. Case 
no. 8 is worse than case no. 1, but case no. 10 is not worse than no. 4, as might 
be expected. 


DYNAMIC INTERACTION STUDY SUMMARY 
[Fundamental frequency in pitch - 0.046 Hz] 





Oscillations 



Maneuver 

Fuel 

Payload 

Vehicle 

ROT-JET-CMD 

comments 

1 

5° Roll 

6.1 

stable 

stable 



2 

Pitch 

4.2 

stable 

stable 



3 

Yaw 

6.3 

stable 

stable 



4 

5° Roll 
with att. hold 
RL-.Ol, DB=.l 

12.0 

diverge 

diverge 

roll , pitch 

results after 
clampdown 

5 

Astronaut 
Forced Osc. 

- 

stable 

stable 


results after 
autopilot on 

6 

Jet F3U 
Stuck on 

- 

stable 

stable 



7 

Jet R4U 
Stuck on 

j 

stable 

stable 

yaw 


8 

3.2° Roll 
Rate =.2 

' 6.1 

stable 

stable 

roll up to 
60 sec 

.046 Hz pulse 

9 

5° Pitch 
with att. hold 

7.5 

stable 

stable 

pitch 


10 

3.2° Roll 

with att. hold 

10.5 

stable 

stable 

roll until 
clampdown 

i 
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This figure and the next one are examples of the output from two SLS runs. 
Note that in both cases oscillatory commands are being sent to the jets, but only 
in the roll case are diverging oscillations seen. 



5 deg pitch iriih rl = deg/iec, db s 1 deg 
with an attitude hold at 60 sec with 
rl “ .OX dcg/sec, db « 1 deg 
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6 dfig roll with rl » .02 deg/sea db ■« 1 deg 
wHh an aiiiiude hold ei 60 vec with 
rl e .01 deg/eec, db = 1 deg 
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The Three Degree of Freedom Payload SImulatlan 
for Orblter/Payload Parametric Studies 


The General Dynamics Beam 


• Base Properties 


• payload/Ij^j^ vehicle = .698 

• ^roll 

.06^12 Hz 

^Pltch 

.049 Hz 

^yaw 

.667 Hz 


• Results with base model 

• payload diverges In pitch 

• oscillatory Jet firings after rate limit change 

• Trends with parameter changes 

• vehicle oscillations are greatest at ff-Qij = .0^16 Hz 

• behavior Insensitive to yaw frequency and Payload changes 

• Increasing rate limits and deadbands can limit small 
oscillations 

• vehicle oscillations Increase as payload inertia Increases 


The last four figures give results from the parametric study. The first 
case is for a three-degree-of-f reedom model of the GD beam. The second is for 
a payload with a roll bending frequency of 0.02 Hz and a larger inertia ratio. 
In the latter case a more severe oscillatory rotation command is seen and the 
vehicle maneuvering is more sluggish. 
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5® Roll Maneuver 
rl - .02 ®/sec . db - .5® 
ot 60 sec set 
rl ■ .01 ®/sec ^ db - .1® 

with 

General Dynamics Beam as the Payload 
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RJC Roll RJC Pitch RJC Yaw 












LARGE SPACE STRUCTURES CONTROLS RESEARCH AND DEVELOPIIENT 


AT MARSHALL SPACE FLIGHT CENTER - 
STATUS AND FUTURE PLANS 


H. Buchanan 

Marshall Space Flight Center 
Systems Dynamics Laboratory 
Huntsville , Alabama 
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OBJECTIVE I: STABILITY AND MODAL CONTROL 

DEMONSTRATE THAT THE FIRST NINE MODES (THREE RIGID + SIX FLEX) 

OF THE SEPS TEST ARTICLE CAN BE ACTIVELY CONTROLLED. 

FEATURES: 

• LOW FREQUENCY (f < 1 Hz). 

• ACTIVE MODAL DAMPING - EXPERIMENT GOAL OF 10%. 

• CONTROL OF ASYMMETRIC STRUCTURE WITH COUPLED MODES. 

• INVESTIGATE EFFECT OF CONTROLLER SATURATION ON DYNAMICS. 


SEPS SOLAR ARRAY FLIGHT TEST MODES 


FREQUENCY Hz 

DESCRIPTION 

0 

RIGID BODY 

0 

RIGID BODY 

0 

RIGID BODY 

.032 

OUT OF PLANE* BENDING 

.035 

IN PLANE BENDING + TORSION** 

.059 

IN PLANE BENDING + TORSION 

.096 

OUT OF PLANE BENDING 

.117 

IN PLANE BENDING + TORSION 

.165 

OUT OF PLANE BENDING 
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•PLANE OF SOLAR BLANKET 
•TORSION OF MAST 
















OBJECTIVE: TO DEVELOP A MULTILEVEL CONTROL APPROACH WHICH SUPPORTS A MODULAR OR BUILDING BLOCK 

APPROACH TO THE BUILDUP OF SPACE PLATFORMS. 



OUTLOOK: CONCEPT HAS BEEN DEVELOPED AND TESTED IN THREE-AXIS COMPUTER SIMULATION UTILIZING A FIVE-BODY 
MODEL OF A BASIC SPACE PLATFORM MODULE. ANALYTICAL EFFORTS HAVE CONTINUED TO FOCUS ON 
EXTENSION OF THE BASIC THEORY AND SUBSEQUENT APPLICATION. 
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DEPLOYABLE ANTENNA SURFACE SHAPE CONTROL 



o OBJECTIVE - DEVELOP PRELIMINARY SPECIFICATIONS 
FOR A FLIGHT EXPERIMENT TO EVALUATE SEVERAL 
ALGORITHMS FOR CONTROLLING THE SHAPE OF LSS. 

o STATEMENT OF WORK SUMMARY 

- DEMONSTRATE ANALYTICALLY THE FEASIBILITY 
FOR SUCH AN EXPERIMENT. 

- SPECIFY HARDWARE AND SOFTWARE REQUIREMENTS. 

- IDENTIFY REQUIREMENTS WHICH WOULD IMPACT 
CURRENT DESIGN. 

- DEFINE A FLIGHT TEST PLAN. 


SAFE CONTROL EXPERIMENT 
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DISTURBANCE ISOLATION AND LOAD ALLEVIATION DURING MANEUVERS 
DEMONSTRATE THAT DISTURBANCES ORIGINATING IN THE 
ORBITER (VCS FIRINGS AND CREW MOTION) CAN BE EFFECTIVELY 
ISOLATED FROM THE TEST ARTICLE BY MEANS OF SOFTWARE AND 
ACTIVE CONTROL. IN A SIMILAR MANNER LOADS IMPOSED ON THE 
STRUCTURE BY MANEUVERING WILL BE ALLEVIATED. 

• 1 X 10~*3 g DISTURBANCE LEVEL. 

• 33 NM TORQUE ALLEVIATION. 

POINTING 

DEMONSTRATE CONTROLLER CAN POINT BASE OF APPENDAGE TO 1 
sec ACCURACY (EXPERIMENT GOAL). NO STAR TRACKER SUN SENSOR 
INPUT WILL BE USED AND PERIOD OF PERFORMANCE WILL BE SHORT TO 
MINIMIZE RATE GYRO DRIFT. 
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Research on the Control of Large Space Structures* 


E. D. Denman 

Department of Electrical Engineering 
University of Houston 
Hou s ton , Te xa s 


♦Paper not presented at conference. T^iis work was partially supported by NASA 
Langley Research Center under grant NSG-1603. 
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INTRODUCTION 


The research effort on the control of large space structures at the University 
of Houston has concentrated on the mathematical theory of finite-element models; 
identification of the mass, damping, and stiffness matrix; assignment of damping to 
structures; and decoupling of structure dynamics. The objective of the work has been 
and will continue to be the development of efficient numerical algorithms for analy- 
sis, control, and identification of large space structures. The major consideration 
in the development of the algorithms has been the large number of equations that must 
be handled by the algorithm as well as sensitivity of the algorithms to numerical 
errors. 

The finite-element model that has been used in the linear second-order matrix 
differential equation 


Md^xC t). + Cdx(tl ^ (1) 

dt^ 


where is the mass^atrix, CeR™^ is the damping matrix, KeR^^’^^^^s the 

stiffness matrix, x(t)eR"' is the node displacement vector, and f(t)eR is the 
forcing function vector. 

The Laplace transform of equation (1) gives the matrix equation 


[Ms^+Cs+K]x(s) = B(s) 


( 2 ) 


where B(s) contains the initial condition information as well as the forcing func- 
tion. If s is replaced in equation (2) by X, equation (2) then takes the form of 
a lambda matrix: 


A(A)X(X) = B(A) 


( 3 ) 


If it is assumed that the initial conditions are zero and no forcing function is 
present, then 


A(X)X(X) = [MX^+CX+K]X(X) 


( 4 ) 
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is the homogeneous equation that will be of interest. The latent roots of 

A(X) are given by 


det A. (A) = det [MX^+CX+K] » 0 


( 5 ) 


and the latent vectors are obtained from 

A(Xpy^ = [MxJ+CX^+K]y^ = 0 


provided that all latent roots are distinct, of multiplicity one. 

Lancaster (ref. 1) and Dennis et al. (ref. 2 ), as well as others, have published 
material on lambda matrices. One of the tasks during the past research period has 
been to consider their work as well as extensions of the algebraic theory of lambda 
matrices for the control of structures. 


ALGEBRAIC THEORY OF LAMBDA MATRICES 

A comprehensive treatment of the algebraic theory of lambda matrices cannot be 
presented in this short paper. Only the essentials necessary to understand the damp- 
ing assignment problem will be given. 

Consider the lambda matrix 

A(\) = Q[IX^+CX+K]Q^ = QA(X)Q^ (6) 


T 

where Q is the Cholesky matrix of the decomposition of M = QQ . Hie normalization 
process in equation (6) will retain the symmetry of C and K. The lambda matrix 
A(\) will be considered in the following work and it will be assumed that the latent 
roots of A(\) appear in conjugate pairs as well as being distinct. This restric- 
tion is not necessary but is made to simplify the analysis. 

It can be shown that the matrix 


A 


c 



(7) 
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has eigenvalues 
eigenvectors of 


that are equal to the 
contain the right and 


latent roots of A ( \) , 
left latent vectors of 


The right and left 
A(X) with 


and 



'ci 


(A^I+C)z^ 



( 8 ) 


( 9 ) 


where and are the right and left eigenvectors, respectively. An eigen- 

projector, P^, will be defined as the matrix 


P 


i 


T 

^ci^ci 


^ci^ci 


( 10 ) 


where the matrix A^^^ defined by 


A 


ci 



A = A 
c Cl 


( 11 ) 


will have the Scume eigenvectors as but all eigenvalues will be zero except for 

which will be the same as in A^. "The eigenprojectors are the matrix residues of 
the partial fraction expansion of [A^(X)]“^ with 


[A^(A)] 


2m 


I 

1=1 



( 12 ) 


where 




(13) 


224 



A latent projector will 
tion expansion of [A(\)]“^ 


be defined as the matrix residues of the partial frac- 
where 


with 


[A(X)]“^ 


2m 


^1 * xi”, 


(14) 


(15) 


The latent projectors are also given by 








dA(X^) 

dX 


(16) 


where and are the latent vectors of A( X) . 

To complete the limited discussion of the algebraic theory, the relationship 

that exists between the eigenprojectors and latent projectors must be given. By 

inverting A (\) it can be shown that 
c 


P 


i 



X^I+C I 
-K 


(17) 


with the past definitions holding. 


ASSIGNMENT OF DAMPING 

An algorithm to assign damping to a particular undamped mode has been developed, 
although the theory has not been fully explored at this time. Considerable work 
remains to be carried out to develop and test a comprehensive algorithm. Therefore 
the material in this section is only preliminary. 

Consider the undamped structure with the associated lambda matrix 

A(X) = IX^+K (18) 
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If it is assumed that K is symmetric and positive definite then all latent roots of 

A(\) lie along the jo) ^xis. Since K is positive definite all of its eigenvalues 

will be real with = w and the eigenprojectors of K are given by 

1 ni 


_ - m P. 

(AI-K)"^ = I 

i=l (X-X^) 


(19) 


where it can be shown that is given by 


P = P 
^i i 




( 20 ) 


which will be <3efined as the augmente<3 latent projector. It is not difficult to show 
that the undamped is 



( 21 ) 


but since the undamped matrix is 


A = 

cu 


0 I 


-K 0 


then 


m 


^ 2 *2 
I ) - 0 

1-1 1 1 


m 


/s. + 

I (PiXi+p^x^) - I 

1=1 


( 22 ) 


(23) 

(24) 
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The damped case has a similar form with the exception that the upper diagonal block 
is modified with 



C25) 


Suppose now that the undamped structure is to have damping added to the 

it 

jth mode (i.e., and \j) and all other modes are to remain undamped. Using 

equations (21) and (22), A becomes 


A 

cd 


m 


m 


r* ^2 *2\ 

I (P/^+P X ) 
i=l 


(P X +P X )K 

1=1 


m 


r r'' , 

(P^X^+P^X^) 

5f /s 2 *2^ 

I (P Xj+P X ) 
1=1 1111 


(P.X>P.X/+CA.+C X ) p.A.+P.X. 

jjjj 33 JJ 33 J3 


/\ * 

.(PjXj«^Xj)K 


which can be rewritten as 


A ,= 

cd 


0 

m 

Y 


L» 

1=1 



m 


-1 ?iK 
1=1 ^ 

0 

Li?^j 



>N 2 *2 

p,a:+p,a/ 
i 3 i i 


-P.K 

j 


-c. 

3 


(26) 


(27) 
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and also in the form 


cd 


= A 


cd 


2m 

I 

i=l 




A ,P. 
cd j 


a'^+a. 


(28) 


The matrix denotes the matrix constructed from the sxmimation and Aj is the 

complement to A"*”. 


thus 


The last step in the algorithm is to recognize that A"*” is orthogonal to Aj ; 


p m - 


” — 


• - 

0 I Pi 

i=l ^ 


0 Pj 


0 0 






m 





-7 P^K 0 


-P.K -C. 


0 0 

L i=i ^ J 


- j J J 


m « 




(29) 


m 


Finally, C- is the matrix that is orthogonal to I eigenprojector 

^ i=l ^ 

Pj must be orthogonal to all as the set of eigenprojectors have the orthogonal 

property 


P P 

i j 


= 0 


i ^ j 


(30) 


The Cj matrices are therefore nothing more than a.P. where a. is a scaling 
constant. The trace of Cj is ^ ^ ^ 


trace (C^) - = 2a^ (31) 
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where trace (P . ) 
of A(X) with ^ 


1. ;itie scaling constants are then used to place the latent roots 


a 


i 



( 32 ) 


2 

i 



( 33 ) 


An example will now be given to illustrate the algorithm, 
select K as 


Let C = 0 and 


K - 


300 -200 

-200 350 


which has eigenvalues ±j 11. 1105 and ±j22.9468. Ttie augmented projectors are 


and 



0.562017 0.496139 

0.496139 0.437983 



0.437983 - 0.496139 

- 0.496139 0.562017 


The matrices 


C 


1 


1 0.8827824 

0.8827824 0.7793048 
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and 


^2 = 


1 -1.1327785 

-1.1327785 1.283187 


are orthogonal to and p^ 

by constants and remain orthogonal to P^ and 

The eigenvalues of 


Note that the columns in and C2 can be scaled 


A = 
c 


-K -C, 


are given in the computed results below, 


A(I,J) MATRIX 


0 

0 

-300 

200 


0 

0 

200 

-350 


1 

0 

0 

0 


0 

1 

0 

0 


EIGENVALUES 

1 

2 

3 

4 


REAL 

0 

0 

3.04403E-22 

3.04403E-22 


IMAGINARY 

22.9468 

-22.9468 

11.1105 

-11.1105 


No damping 


0 

0 

-300 

200 


A(I,J) MATRIX 


0 

0 

200 

-350 


1 

0 

-1 

-.882783 


0 

1 

-.882783 

-.779305 


EIGENVALUES 


REAL 


IMAGINARY 


1 

2 

3 

4 


3.99893E-08 

3.99893E-08 

- .889652 

- .889652 


22.9468 

-22.9468 

11.0748 

-11.0748 


Damping in lowest mode 
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A(I,J) MATRIX 


0 

0 

-300 

200 


0 

0 

200 

-350 


1 0 

0 1 

-1 1.13278 

1.13278 -1.28319 


EIGENVALUES 

1 

2 

3 

4 


REAL 

-1.14159 

-1.14159 

2.84795E-08 

2.84795E-08 


IMAGINARY 

22.9184 

-22.9184 

11.1105 

-11.1105 


Damping in highest mode 


A(I,J) MATRIX 


0 

0 

-300 

200 


0 

0 

200 

-350 


1 0 

0 1 

-1 .125 

.125 -1.03125 


EIGENVALUES REAL 


IMAGINARY 


1 

2 

3 

4 


-.570798 

-.570798 

-.444827 

-.444827 


22.9399 

-22.9399 

11.1019 

-11.1019 


Damping in both modes 


The undamped system matrix and its eigenvalues are given first. The lowest mode 
was then damped with the system matrix and its eigenvalue given. The third test was 
to include damping for the highest mode where C = C2 . The last test was to combine 
the damping for the two modes with 


■ 2 «=l+‘=2> ■ 


c^Pl+ajPj 


0.889655P + 1.141595P 


The system matrix and its eigenvalues are then given. It should be noted that the 
new values of Oj of the two modes are 


0.889652 

Ui 2 


0.444827 
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and 


1.141595 
02 2 


.570797 


The algorithm allows dcunping to be assigned to all modes and to the prediction of the 
location of the system eigenvalues. 

The major problem that remains to be resolved is that of constructing the C 
matrix, which is symmetric and positive definite as well as realizable. It will be 
assumed that C and K are tridiagonal with the off-diagonal elements having a 
smaller magnitude than the diagonal elements. Hie off-diagonal elements must be 
negative with positive diagonal terms. 


IDENTIFICATION OF M, C, AND K 

The quadrature algorithm for identifying the mass, damping, and stiffness 
matrices of a structure is still under study. The algorithm performs well with 
simulated data, but attempts to utilized data from the beam experiment at NASA 
Langley Research Center have not produced useable results. It is believed that the 
conditions on the beam during the data collection may not have satisfied the require- 
ments of the algorithm. 

The quadrature identification, as well as other algorithms, will receive major 
attention in the future. One of the major problems in developing an efficient 
algorithm is the availability of test data, either by computer simulation or from a 
test bed. Hie Langley Research Center experimental beam is a suitable test bed, but 
planning and running test data are costly and time consuming. Plans for the future 
include construction of an electronic analog test bed that will be low cost and will 
provide flexibility in the types of structures that can be simulated. Hie test bed 
wiilll, be used primarily in the development of the identification algorithm, but it can 
also be used for developing control algorithms. 


SIMULATION OF STRUCTURES 

Some of the preliminary work on designing a structures test bed at the 
University of Houston has been completed. Hie type of simulator considered is an 
electronic analog that will be constructed from low-cost operational amplifiers, 
resistors, and capacitors. Data acquisition will be handled by a PDP 11/70 computer 
for signal processing. Hie PDP 11/70 is equipped with analog-to-digital and digital- 
to-analog converters and has adequate storage to perform on-line data collection and 
processing. 


CONCLUSIONS 

Some of the work that has been supported under grant from NASA Langley Research 
Center has been described in this report. Preliminary details on the assignment of 
damping are covered, as well as some information on identification and simulation. 
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The identification work has not been described since the algorithm is given in refer- 
ence 3. The simulation facility is still in the planning stage and a decision to 
build such a test bed has not been made. Further descriptions of this research are 
given in references 4 and 5, 
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ACTIVE CONTROL OF A FLEXIBLE BEAM 


Because of Inherent low damping and high flexibility, large space structures 
may require some form of active control of their dynamics. Because of the apparent 
inability to accurately model the dynamics of these structures, methods for 
parameter adaptive control are now being developed at Langley, The process 
currently being studied is shown in the block diagram below. This approach uses a 
digital computer to process discrete sensor data. Identify modal parameters, 
calculate modal control gains, and then convert the modal forces to real forces. 

The last two blocks are the topic of this presentation. Some of the problems con- 
sidered are: Cl) the possibility that there loay be many modes to control with 

limited amounts of hardware, and (2) the reqxilred accuracy of identified structural 
parameters. 


BACKGROUND: 

- NEED TO CONTROL FLEXIBLE MOTION OF LARGE SPACE STRUCTURES 

- ABILITY TO ACCURATELY MODEL THE DYNAMICS OF THESE STRUCTURES IS UNCERTAIN 

- THEORY NOW BEING DEVELOPED FOR PARAMETER ADAPTIVE CONTROL OF THESE STRUCTURES 



PROBLEMS: 

- POSSIBLY MANY MODES TO CONTROL 

- LIMITED HARDWARE (COMPUTATION. SENSORS. ACTUATORS) 

- HOW WELL STRUCTURAL PARAMETERS MUST BE IDENTIFIED FOR CONTROL 
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BESEARCH TASK 


The specific research task was to design a digital control scheme to suppress 
vibration of a homogeneous free-free beam. A digital computer simulation algorithm 
was then used to test (1) the effects of controlling more modes than available 
actuators, and (2) the sensitivity to identified structural parameters. 


• DESIGN A DIGITAL CONTROL SCHEME TO SUPPRESS VIBRATION OF A 
HOMOGENEOUS FREE-FREE BEAM 

• EXAMINE EFFECT ON STABILITY OF: 

- FEWER ACTUATORS THAN CONTROLLED MODES 

- ERRORS IN STRUCTURAL MODEL PARAMETERS 


• TEST WITH AN EXISTING SIMULATION ALGORITHM 



MATHEMATICAL MODEL 


The continuous beam was modeled by using the SPAR finite element a^lgorlthm 
which generates mode shapes and frequencies. These were used to write a modal 
representation of the beam dynamics which was used to design the control gains. 






CONTINUOUS nODEL 



FINITE ELEMENT MODEL OF BEAM 

- MODE SHAPES AND FREQUENCIES FOR 25 ELEMENT MODEL ARE 
GENERATED BY SPAR 

[M]£ + (K]x * F 



• MODAL REPRESENTATION 


X = [E] q 

J q + J a “ U J li 

- SET OF 50 UNCOUPLED 2ND-0RDER SYSTEMS 


[E]'F 


• CONTROL SYSTEM DESIGNED USING THIS MODAL REPRESENTATION OF THE 
STRUCTURAL DYNAMICS 
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DISCRETE TIME MODEL 


In order to simplify the digital simulation and prepare for eventual digital 
implementation, the modal equations of motion were discretized. This results in 
the scalar equation which shows the present modal anqplitude to be a function of the 
past two amplitudes and the past two controls. There is a discrete time transforma- 
tion analogous to the Laplace transform which results in a characteristic polynomial 
in z. The roots of this polynomial may be plotted in the con^ilex plane with 
stability represented by magnitudes of less than 1. Analysis of the control system 
will be done primarily in this z-plane. 


• NECESSARY FOR: 

- DESIGN OF DIGITAL SIMULATION 

- EVENTUAL DIGITAL IMPLEMENTATION 

• DISCRETE EQUATION OF MOTION: 

''(k) " ^1 '’(k-1) ^ ^2 '’(k-2) ^ ^1 ‘^(k-1) “(k-2)' 

where A = 1‘(w ,t), t = sampling interval, 

B = f(w ,e,x), £ = mode shapes. 



I »—■ ' 

(k- 1 )t 


• DISCRETE TIME TRANSFORM (ANALOGOUS TO LAPLACE TRANS.) RESULTS 
IN A CHARACTERISTIC POLYNOMIAL IN Z: n . 


S-PLANE 



-PLANE 


' i 

^ UNSTABLE 
REGION 


UNSTABLE 

REGION 


CONTINUOUS TIME 


DISCRETE 


Time 
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CONTBOL DESIGN APPEOACH 


The nodal control design approach Is to choose desired closed-loop roots from 
which the modal controller gains can be calcxilated. The modal control forces nay 
be calculated directly, and the actual control forces can be calculated using a 
pseudo-inverse . 


• CHOOSE DESIRED CLOSED LOOP ROOTS FOR EACH MODE 

• CALCULATE MODAL CONTROLLER GAINS 

• CALCULATE MODAL CONTROL FORCES 

• CONVERT MODAL FORCES TO ACTUATOR FORCES 
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DIGITAL CONTROL OF ONE MODE 


The control of one mode is achieved by using the minimum order control law 
required for pole placement. This is of the same form as the modal equation of 
motion. The closed loop controller has a fourth-order characteristic equation as 
shown. The coefficients of this equation are determined from the desired closed- 
loop roots and are functions of the mode and control coefficients in the plant and 
control modal equations. The control objective is to achieve the desired roots 
defined by (a^b^Cyd} by solving for the controller gains C2» D^, ^ 2 ^* 


• MINIMUM ORDER CONTROL LAW F0R^P0LE PLACEMENT 

“(k) " ^1 ‘’(k-1) * ^2 ‘’(k-2) * °1 “(k-1) ^2 “(k-2) 

- NOTE SAME FORM AS PLANT MODEL 

f‘*(k) " ^1 ‘’(k-1) ^2 ‘'(k-2) ‘^(k-1) ®2 “(k-2)^ 

- C'S & D'S ARE FOUR CONTROL GAINS 

• CLOSED LOOP CONTROLLER 



- HAS DISCRETE TIME CHARACTERISTIC EQUATION 

z^+az^+bz^+cz+d =0 

where {a ,b ,c ,d } » f(A ,B ,C ,D ). 

• CONTROL OBJECTIVE: 

- ACHIEVE DESIRED CLOSED LOOP ROOTS (AS DEFINED BY a.b.c.d) 
BY CALCULATING THE CONTROL GAINS. 
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CLOSED-LOOP CHARACTERISTICS 

The control problem is shown In the block diagram below. The parameter 
uncertainties affect the calcvilatlon of the modai controller gains and a pseudo- 
inverse results in a least squares type error in the actual forces applied to the 
beam. A typical set of closed-loop roots is shown in each of the two complex 
plane plots. The design criterion is to place the roots of the controlled modes on 
a constant damping line in the s-plane. This line maps onto a constant-radius 
circle on the cos^lex z-plane, with uncontrolled roots on the unit circle. 


0 CONTROL PROBLEM 



0 REPRESENTATION OF CLOSED LOOP DYNAMICS 

Z-PLFINE 
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RESULTS - CASE 1 


A baseline set of results is shown below. Exact parameters are used to 
calculate the control ^ains and five modes are controlled with five actuators. 
Note that all roots are calculated exactly so that all modal aiq> litudes have the 
same decay envelope. Also, mode six, which is not controlled, does show "minor” 
excitation and continues to "ring" after control to the other modes is stopped. 
This is not evident from the figure. 


• CASE I : EQUAL NUMBER OF ACTUATORS AND CONTROLLED MODES 


S-PLRNE 



TlfC. Stc. 5CC. 
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BESULTS - CASE II 


Thla next case illustrates an atteai5)t to control anore modes than available 
actuators. Here it is no longer possible to solve exactly for the desired roots, 
and one root does become nnstable, as shown by the increasing anplitude of the 
fo^lrth mode. 


• CASE .1 1 : FEWER ACTUATORS THAN CONTROLLED MODES - UNSTABLE 


5 CONTROLLED 


Z-PLRNE 
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RESULTS - CASE III 


In the third case it Is shown that attempting to control more modes than 
available actuators does not necessarily mean the system will be tinstable. Note 
here that the calculated roots are closer to the desired roots and that, while the 
fourth mode is near the unit circle, it now is slightly stable. This can be 
explained by considering that the additional actuator and mode provide one more 
data point for the least squares fit. 


f CASE III: FEWER ACTUATORS THAN CONTROLLED MODES - STABLE 


6 

5 


CONTROLLED 

ACTUATORS 


Z-PLRNE 


IH 



mode H n 

amplitude ° 


-I 


y\/VV WW\MA/'u'VVWA^w 
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PARAMETER ERROR SENSITIVITY 


The results of the parameter error sensitivity study are shown below. Errors 
was placed on the modal frequency and damping parameters of the modes selected for 
control. Experimental results show that the parameter estimator may have errors on 
the order of ±5 percent. The locus of roots calculated using a parameter error 
range of ±20 percent is plotted below. The magnitude of the error at which the 
modes became unstable is summarized in the table. Note that as the mode number 
increases, the sensitivity decreases for the first three modes. The high sensi- 
tivity of mode six in case III is unexplained at this time. 


f 


ERROR ON 

NO ERROR 
BASED ON 


MODAL FREQUENCY AND DAMPING PARAMETERS ONLY 

^ “(k-i) ^ 

ON RIGID BODY PARAMETERS 

EXPERIMENTAL PARAMETER ESTIMATE ERRORS OF + 5 % 


5 CONTMUED 
5 ACTUATORS 

Z-PLRNE 




jPARAMETER ERROR AT INSTABILITY 

CASE 

'mode 3 

M0r>E 4 { M03E 5 

MODE 6 

1 5 CONT. 

-10 % 

-12% 

-14% 


5 ACT. 

+8 % 

— 



I I 5 CONT. 

-6 % 

0.0 % 

-14 % 


R ACT. 

+ 14% 

— 

— 


Ill 6 CONT. 

-8 % 

-12 % 

-16% 

-2 % 

5 ACT. 

+8 % 

— 

— 

+2 % 
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CONCLUSIONS 


1. LIMITED ACTUATORS: 

- THIS DESIGN PROCESS YIELDS UNDESIRABLE CLOSED LOOP DYNAMICS 
WHEN THE NUMBER OF CONTROLLED MODES EXCEEDS THE NUMBER OF 
AVAILABLE ACTUATORS. 

2. PARAMETER ERROR: 

- ERRORS WITHIN THE RANGE OF EXPERIMENTAL RESULTS CAN CAUSE 
INSTABILITY. 

- CONTROL SYSTEM MUST BE MADE MORE TOLERANT OF PARAMETER ERROR. 
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OPTIMAL DAMPING FOR A 


TWO-DIMENSIONAL STRUCTURE 


W. D. Pilkey and B. P. Wang 
University of Virginia 
Charlottesville , Virginia 



SYNTHESIS OF THE DAMPING MATRIX 
FOR SPECIFIED DAMPING RATIOS 


[M]{x) + tkHx) -(F) {X) is n X 1 

CONTROL FORCES 


(F) - 

[B]{u) 

{u} 

is 

A 

X 1 



A “ number 

of 

actuators 

{U} = 

-Id] {X) 

IDJ 

is 

A 

X n 

or 






{f} = 

= -tc]{x} 






where 

[C] = [B][D] 
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PROBLEM 


Find such that 


J - X: » C l 

irj ^ 


is minimized subject to the constraints 


C >C for£“ltoL 

t ip 


and 


<=u > « ■ ‘=13 - ‘=31 


» prescribed damping ratios for the ith mode 


Eigenvalue Problem: 


(3 [M] + s [C3 + (K])tP> - (O) 


with roots 


s , 


'A - 


i in 


in 
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DAMPED EIGENVALUE PROBLEM-. 



ALTERNATE EFFICIENT 
FORMULATION FOR 
CHARACTERISTIC EQUATION •• 

dei([I] + s[R(s)][^]) s 0 

ASSUMING NO. OF DAMPERS « 
NO. OF 0.0. E 


[R] u SUBMATRIX OF [R] 


SPECTRAL REPRESENTATION 



G|(s‘■^uf) 


y 011 <Pm 



DERIVATION OF CHARACTERISTIC 
EQUATION : 



CHARACTERISTIC EQ. IS ; 

<Tet([I]-^s[^][e])=0 


SPECIAL CASE .■ 

INTRODUCE SINGLE DAMPER 
AT D.O.F. J 

CHARACTERISTIC EQUATION; 

/ + scRjj s 0 

OR 

■L - -I ^ 
c ut 
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DEFINITION : 

OPTIMAL DAMPER LOCATIONS 
FOR A PARTICULAR MODE 
ARC WHERE EITHER 

MAXIMUM DAMPm CAN 
BE INTRODUCED 

OR 

ACHIEVE GIVEN 0AMPIN6 
WITH MINIMUM DAMPING 
CONSTANTS 



TYPICAL ROOT LOCUS PLOT 



MINIMUM CONSTRAINED 
FREQUENCY CRITERION 
(MCFC) 


THE OPTIMAL HAMPER 
LOCATION IS WHERE THE 
CONSTRAINEO FREQUENCY 
IS A MINIMUM 



TYPICAL ROOT LOCUS 
FOR INTRODUCING MORE 
DAMPERS THAN NO. OF 
RIGID BODY MODES 



MAXIMUM FREQUENCY 
SEPARATION CRITERION 
(MFSC) 


THE OPTIMAL HAMPER 
LOCATION IS WHERE 
THE SEPARATION BETWEEN 
THE CONSTRAINED 
FREQUENCY AND THE 
NATURAL FREQUENCY 
IS MAXIMUM 


SVNTHESIS QE SINSIE DAMPER 
CHARACTERISTIC EQUATION IS 
f(s) = / + sRjf(s)c 

• / + cPIs) + JcQtt) mO 

SPECIFY ( (OR CJ^) 

THEN , 

S, = 

SOLVE FOR (OR^) 

FROM Qls,)uO 

THE DESIRED DAMPIN6 CONSTANT 

c s-1. 

P<H) 

C>0 IF s, IS REAUZABLE 



SVNWSSIS OF MULTIPLE 
DAMPERS 


CHARACTERISTIC EQUATION IS 
fis) X dti (I + sR(s) S) 

- P(s,£)+JQ(s,£)xO 

CASE I 

SPECIFY nc $ 'S 
THEN SOLVE FOR 


£ -[ Cqf 

AND 

> ^0^ >*"J ^/lie 

FROM THE Znc EQS 

P(h,£) ’0 


Q(Si,£)xO 


FOR (•! TOnc 


CASE 2 

SPECIFY E wm E<nc 

THEN C AND 

CAN BE FOUND FROM 

MINIMIZE JsZCi 

/«/ /«/ 


mTH q:»0 

«. 9 O 



RANK OF SOME TWO PAMPER LOCATIONS 
ACCORDING TO NiCFC 


OAMFER 

LOCATIOMS 

CMIteSPMOINO 
LOWEST COHSTEAmtO 
FREOUEMCT 

RANK 

/.«< 

4S 8 RAOfSec 

/ 

I.SO 

e\B9 

3 

II, SB 

7. so 

4 

34,44 

6.S4 

Z 


DAMPING CONSTANTS REQUIRED TO 


ACHIEVE e^:^0.6 


\OAMKR LOCATION 

c, 

REQUlAiO 

Ct c,*c^ 

NO. 1 

NO. 2 

1 

86 

a 296 

o.m 

0.432 

1 

SO 

0.262 

1.024 

1.266 

II 

50 

0.311 

1.130 

1.441 

34 

44 

0.496 

0.565 

1.064 


THREE DAMPER EXAMPLE 


DAMPER 

LOCATIOMS 



Os 

\Us-Os\ 

1, 50, ea 

t, 44, 78 

0,62 

6.40 

13.13 

7.35 

11.56 

8.48 

4.81 

7.83 


MAXIMUM 0AMPIN6 ACHIEVeO IN FIRST 

TWO VieRATORT MOORS 

{WITH SAME OAMFINS COHSTANT AT ALL 

LOCATIONS) 


SANDER 

LOCATIONS 

(4 

OAMPBR 

VALUE 


OAMPEK 

VALUE 

1,50,88 
!, 44, 78 

0.95 

0.60 

0.34 

0.31 

0.05 

0.46 

0.48 

0.34 




CONCLUSIONS 


cm Tern for sflection of 

OPTIMAL DAMPER LOCATIONS 
PRESENTED 

DAMPING SYNTHESIS PROBLEM 
FORMULATED AND APPLIED 
TO NASA GRILLAGE MODEL 


DAMPING SYNTHESIS 
EXAMPLE I 

DESIRED ^ s a£ 

= O.S 


DAMPER LOCATIONS 


OAF. 1,44 

,7a 

RESULTS : 

LOCATION 

GAIN 

/ 

0.272 

44 

0514 

78 

0.269 










DAMPING SYNTHESIS EXAMPLE 2 


DESIPED = 0.7 
= 0.S 

USE G DAMPERS AT D.O.F . 
1,11,39,50,78^86 


RESULTS : 

LOCATION! Ci Cl 


i 

0,Z! 

0.240 

It 

0.25 

0.240 

39 

0.20 

0. 121 

SO 

0.20 

0.131 

78 

0.20 

0.247 

88 

0.21 

0.24S 

ACHIEVED 

DAMPING 



(5 * 0.71 

<- 
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